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TECHNICAL FIELD 

The invention relates to flowmeters. 



BACKGROUND 

15 Flowmeters provide information about materials being transferred through a 

conduit. For example, mass flowmeters provide a direct indication of the mass of material 
being transferred through a conduit. Similarly, density flowmeters, or densitometers, provide 
an indication of the density of material flowing through a conduit. Mass flowmeters also 
may provide an indication of the density of the material. 

20 Coriolis-type mass flowmeters are based on the well-known Coriolis effect, in 

which material flowing through a rotating conduit becomes a radially traveling mass that is 
affected by a Coriolis force and therefore experiences an acceleration. Many Coriolis-type 
mass flowmeters induce a Coriolis force by sinusoidally oscillating a conduit about a pivot 
axis orthogonal to the length of the conduit. In such mass flowmeters, the Coriolis reaction 

25 force experienced by the traveling fluid mass is transferred to the conduit itself and is 

manifested as a deflection or offset of the conduit in the direction of the Coriolis force vector 
in the plane of rotation. 

Energy is supplied to the conduit by a driving mechanism that applies a periodic 
force to oscillate the conduit. One type of driving mechanism is an electromechanical driver 

30 that imparts a force proportional to an applied voltage. In an oscillating flowmeter, the 
applied voltage is periodic, and is generally sinusoidal. The period of the input voltage is 



1 



Attorney Docket No.: 02052-079004 



chosen so that the motion of the conduit matches a resonant mode of vibration of the conduit. 
This reduces the energy needed to sustain oscillation. An oscillating flowmeter may use a 
feedback loop in which a sensor signal that carries instantaneous frequency and phase 
information related to oscillation of the conduit is amplified and fed back to the conduit using 
5 the electromechanical driver. 

SUMMARY 

The invention provides a digital flowmeter, such as a digital mass flowmeter, that 
uses a control and measurement system to control oscillation of the conduit and to generate 

10 mass flow and density measurements. Sensors connected to the conduit supply signals to the 
control and measurement system. The control and measurement system processes the signals 
to produce a measurement of mass flow and uses digital signal processing to generate a 
signal for driving the conduit. The drive signal then is converted to a force that induces 
oscillation of the conduit. 

15 The digital mass flowmeter provides a number of advantages over traditional, 

analog approaches. From a control perspective, use of digital processing techniques permits 
the application of precise, sophisticated control algorithms that, relative to traditional analog 
approaches, provide greater responsiveness, accuracy and adaptability. 

The digital control system also permits the use of negative gain in controlling 

20 oscillation of the conduit. Thus, drive signals that are 180° out of phase with conduit 

oscillation may be used to reduce the amplitude of oscillation. The practical implications of 
this are important, particularly in high and variable damping situations where a sudden drop 
in damping can cause an undesirable increase in the amplitude of oscillation. One example 
of a variable damping situation is when aeration occurs in the material flowing through the 

25 conduit. 

The ability to provide negative feedback also is important when the amplitude of 
oscillation is controlled to a fixed setpoint that can be changed under user control. With 
negative feedback, reductions in the oscillation setpoint can be implemented as quickly as 
increases in the setpoint. By contrast, an analog meter that relies solely on positive feedback 
30 must set the gain to zero and wait for system damping to reduce the amplitude to the reduced 
setpoint. 
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From a measurement perspective, the digital mass flowmeter can provide high 
information bandwidth. For example, a digital measurement system may use analog-to- 
digital converters operating at eighteen bits of precision and sampling rates of 55 kHz. The 
digital measurement system also may use sophisticated algorithms to filter and process the 

5 data, and may do so starting with the raw data from the sensors and continuing to the final 
measurement data. This permits extremely high precision, such as, for example, phase 
precision to five nanoseconds per cycle. Digital processing starting with the raw sensor data 
also allows for extensions to existing measurement techniques to improve performance in 
non-ideal situations, such as by detecting and compensating for time- varying amplitude, 

10 frequency, and zero offset. 

The control and measurement improvements interact to provide further 
improvements. For example, control of oscillation amplitude is dependent upon the quality 
of amplitude measurement. Under normal conditions, the digital mass flowmeter may 
maintain oscillation to within twenty parts per million of the desired setpoint. Similarly, 

15 improved control has a positive benefit on measurement. Increasing the stability of 
oscillation will improve measurement quality even for meters that do not require a fixed 
amplitude of oscillation (i.e., a fixed setpoint). For example, with improved stability, 
assumptions used for the measurement calculations are likely to be valid over a wider range 
of conditions. 

20 The digital mass flowmeter also permits the integration of entirely new 

functionality (e.g., diagnostics) with the measurement and control processes. For example, 
algorithms for detecting the presence of process aeration can be implemented with 
compensatory action occurring for both measurement and control if aeration is detected. 

Other advantages of the digital mass flowmeter result from the limited amount of 

25 hardware employed, which makes the meter simple to construct, debug, and repair in 
production and in the field. Quick repairs in the field for improved performance and to 
compensate for wear of the mechanical components (e.g., loops, flanges, sensors and drivers) 
are possible because the meter uses standardized hardware components that may be replaced 
with little difficulty, and because software modifications may be made with relative ease. In 

30 addition, integration of diagnostics, measurement, and control is simplified by the simplicity 
of the hardware and the level of functionality implemented in software. New functionality, 
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such as low power components or components with improved performance, can be integrated 
without a major redesign of the overall control system. 

In one general aspect, a digital flowmeter includes a vibratable conduit, a driver 
connected to the conduit and operable to impart motion to the conduit, and a sensor 
5 connected to the conduit and operable to sense the motion of the conduit. A control and 
measurement system connected to the driver and the sensor includes circuitry that receives a 
sensor signal from the sensor, generates a drive signal based on the sensor signal using digital 
signal processing, supplies the drive signal to the driver, and generates a measurement of a 
property of material flowing through the conduit based on the signal from the sensor. 

10 Embodiments may include one or more of the following features. The meter may 

include a second sensor connected to the conduit and operable to sense the motion of the 
conduit. In this case, the control and measurement system is connected to the second sensor 
and receives a second sensor signal from the second sensor, generates the drive signal based 
on the first and second sensor signals, and generates the measurement of the property of 

15 material flowing through the conduit based on the first and second sensor signals. The 
control and measurement system may digitally combine the first and second sensor signals 
and generate the drive signal based on the combination of the sensor signals. 

The control and measurement system may generate different drive signals for the 
two drivers. The drive signals may have, for example, different frequencies or amplitudes. 

20 The digital flowmeter also may include circuitry for measuring current supplied to 

the driver. The circuitry may include a resistor in series with the driver and an analog-to- 
digital converter in parallel with the resistor and configured to measure a voltage across the 
resistor, to convert the measured voltage to a digital value, and to supply the digital value to 
the control and measurement system. 

25 The digital flowmeter also may include a first pressure sensor connected to 

measure a first pressure at an inlet to the conduit and a second pressure sensor connected to 
measure a second pressure at an outlet of the conduit. Analog-to-digital converters may be 
connected and configured to convert signals produced by the first pressure sensor and the 
second pressure sensor to digital values and to supply the digital values to the control and 

30 measurement system. Temperature sensors may be connected to measure temperatures at the 
inlet and outlet of the conduit. 
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The control and measurement system may generate the measurement of the 
property by estimating a frequency of the first sensor signal, calculating a phase difference 
using the first sensor signal, and generating the measurement using the calculated phase 
difference. The control and measurement system may compensate for amplitude differences 
5 in the sensor signals by adjusting the amplitude of one of the sensor signals. For example, 
the control and measurement system may multiply the amplitude of one of the sensor signals 
by a ratio of the amplitudes of the sensor signals. 

When the sensor signal is generally periodic, the control and measurement system 
may process the sensor signal in sets. Each set may include data for a complete cycle of the 
10 periodic sensor signal, and consecutive sets may include data for overlapping cycles of the 
periodic sensor signal. The control and measurement system may estimate an end point of a 
cycle using a frequency of a previous cycle. 

The control and measurement system may analyze data for a cycle to determine 
whether the cycle merits further processing. For example, the system may determine that a 
15 cycle does not merit further processing when data for the cycle does not conform to expected 
behavior for the data, where the expected behavior may be based on one or more parameters 
of a previous cycle. In one implementation, the system determines that a cycle does not 
merit further processing when a frequency for the cycle differs from a frequency for the 
previous cycle by more than a threshold amount. The system may determine whether the 
20 frequencies differ by comparing values at certain points in the cycle to values that would 
occur if the frequency for the cycle equaled the frequency for the previous cycle. 

The control and measurement system may determine a frequency of the sensor 
signal by detecting zero-crossings of the sensor signal and counting samples between zero 
crossings. The system also may determine a frequency of the sensor signal using an iterative 
25 curve fitting technique. 

The control and measurement system may determine an amplitude of the sensor 
signal using Fourier analysis, and may use the determined amplitude in generating the drive 
signal. 

The control and measurement system may determine a phase offset for each 
30 sensor signal and may determine the phase difference by comparing the phase offsets. The 
system also may determine the phase difference using Fourier analysis. The control and 
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measurement system may determine a frequency, amplitude and phase offset for each sensor 
signal, and may scale the phase offsets to an average of the frequencies of the sensor signals. 
The control and measurement system may calculate the phase difference using multiple 
approaches and may select a result of one of the approaches as the calculated phase 
5 difference. 

The control and measurement system may combine the sensor signals to produce 
a combined signal and may generate the drive signal based on the combined signal. For 
example, the control and measurement system may sum the sensor signals to produce the 
combined signal and may generate the drive signal by applying a gain to the combined 
10 signal. 

In general, the control and measurement system may initiate motion of the 
conduit by using a first mode of signal generation to generate the drive signal, and may 
sustain motion of the conduit using a second mode of signal generation to generate the drive 
signal. The first mode of signal generation may be synthesis of a periodic signal having a 

15 desired property, such as a desired initial frequency of conduit vibration, and the second 
mode of signal generation may use a feedback loop including the sensor signal. 

In other instances, the first mode of signal generation may include use of a 
feedback loop including the sensor signal and the second mode of signal generation may 
include synthesis of a periodic signal having a desired property. For example, the control 

20 and measurement system may generate the drive signal by applying a large gain to the 

combined signal to initiate motion of the conduit and generating a periodic signal having a 
phase and frequency based on a phase and frequency of a sensor signal as the drive signal 
after motion has been initiated. The desired property of the synthesized signal may be a 
frequency and a phase corresponding to a frequency and a phase of the sensor signal. 

25 The control and measurement system generates an adaptable, periodic drive 

signal. For example, the meter may include positive and negative direct current sources 
connected between the control and measurement system and the driver, and the control and 
measurement system may generate the drive signal by switching the current sources on and 
off at intervals having a phase and frequency based on the sensor signal. The control and 

30 measurement system may generate the drive signal by synthesizing a sine wave having a 
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property corresponding to a property of the sensor signal, such as a phase and a frequency 
corresponding to a phase and a frequency of the sensor signal. 

The control and measurement system may digitally generate a gain for use in 
generating the drive signal based on one or more properties of the sensor signal. For 
5 example, the control and measurement system may digitally generate the gain based on an 
amplitude of the sensor signal. 

The driver may be operable to impart an oscillating motion to the conduit. The 
control and measurement system also may digitally implement a PI control algorithm to 
regulate the amplitude of conduit oscillation. The control and measurement system also may 
10 digitally generate the drive signal based on the sensor signal so as to maintain an amplitude 
of oscillation of the conduit at a user-controlled value. In support of this, the control and 
measurement system may generate a negative drive signal that causes the driver to resist 
motion of the conduit when the amplitude of oscillation exceeds the user-controlled value 
and a positive drive signal that causes the driver to impart motion to the conduit when the 
1 5 amplitude of oscillation is less than the user-controlled value. 

The control and measurement system may include a controller that generates a 
gain signal based on the sensor signal and a multiplying digital-to-analog converter 
connected to the controller to receive the gain signal and generate the drive signal based on 
the gain signal. 

20 When the digital flowmeter includes a second sensor connected to the conduit and 

operable to sense the motion of the conduit, the control and measurement system may 
include a controller that generates the measurement, a first analog-to-digital converter 
connected between the first sensor and the controller to provide a first digital sensor signal to 
the controller, and a second analog-to-digital converter connected between the second sensor 

25 and the controller to provide a second digital sensor signal to the controller. The controller 
may combine the digital sensor signals to produce a combined signal and to generate a gain 
signal based on the first and second digital sensor signals. The control and measurement 
system also may include a multiplying digital-to-analog converter connected to receive the 
combined signal and the gain signal from the controller to generate the drive signal as a 

30 product of the combined signal and the gain signal. 
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The control and measurement system may selectively apply a negative gain to the 
sensor signal to reduce motion of the conduit. 

The control and measurement system also may compensate for zero offset in the 
sensor signal. The zero offset may include a component attributable to gain variation and a 
5 component attributable to gain nonlinearity, and the control and measurement system may 
separately compensate for the two components. The control and measurement system may 
compensate for zero offset by generating one or more correction factors and modifying the 
sensor signal using the correction factors. 

The control and measurement system may calculate phase offsets for the first and 
10 second sensor signals. The phase offset may be defined as a difference between a zero- 
crossing point of a sensor signal and a point of zero phase for a component of the sensor 
signal corresponding to a fundamental frequency of the sensor signal. The control and 
measurement system may combine the calculated phase offsets to produce a phase 
difference. 

15 The control and measurement system may generate the measurement of the 

property by estimating a frequency of the first sensor signal, estimating a frequency of the 
second sensor signal, with the frequency of the second sensor signal being different from the 
frequency of the first sensor signal, and calculating a phase difference between the sensor 
signals using the estimated frequencies. 

20 When the sensor is a velocity sensor, the control and measurement system may 

estimate a frequency, amplitude, and phase of the sensor signal, and may correct the 
estimated frequency, amplitude, and phase to account for performance differences between a 
velocity sensor and an absolute position sensor. Instead of controlling the apparent 
amplitude of oscillation (i.e., the velocity of oscillation when the sensor is a velocity sensor), 

25 the system may control the true amplitude by dividing the sensor signal by the estimated 
frequency. This correction should provide improved amplitude control and noise reduction. 

The control and measurement system may estimate a first parameter of the sensor 
signal, determine a rate of change of a second parameter, and correct the estimated first 
parameter based on the determined rate of change. For example, the system may correct an 

30 estimated frequency, amplitude, or phase of the sensor signal based on a determined rate of 
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change of the frequency or amplitude of oscillation of the conduit. The system may perform 
separate corrections for each sensor signal. 

The digital flowmeter may be a mass flowmeter and the property of material 
flowing through the conduit may be a mass flow rate. The digital flowmeter also may be a 
5 densitometer and the property of material flowing through the conduit may be a density of 
the material. 

The control and measurement system may account for effects of aeration in the 
conduit by determining an initial mass flow rate, determining an apparent density of material 
flowing through the conduit, comparing the apparent density to a known density of the 

10 material to determine a density difference, and adjusting the initial mass flow rate based on 
the density difference to produce an adjusted mass flow rate. The system may further 
account for effects of aeration in the conduit by adjusting the adjusted mass flow rate to 
account for effects of damping. To further account for effects of aeration in the conduit, the 
system may adjust the adjusted mass flow rate based on differences between amplitudes of 

15 the first and second sensor signals. 

The vibratable conduit may include two parallel planar loops. The sensor and 
driver may be connected between the loops. 

The meter may include a power circuit that receives power on only a single pair 
of wires. The power circuit provides power to the digital control and measurement system 

20 and to the driver, and the digital control and measurement system is operable to transmit the 
measurement of the property of material flowing through the conduit on the single pair of 
wires. The power circuit may include a constant output circuit that provides power to the 
digital control and measurement system and drive capacitor that is charged by excess power 
from the two wires. The digital control and measurement system may discharge the drive 

25 capacitor to power the driver, and may monitor a charge level of the drive capacitor and 
discharge the drive capacitor after a charge level of the capacitor reaches a threshold level. 
The digital control and measurement system also may discharge the drive capacitor 
periodically and to perform bidirectional communications on the pair of wires. 

The control and measurement system may collect a first data set for a period of 

30 the periodic signal and process the first data set to generate the drive signal and the 
measurement. The system may collect a second data set for a subsequent period of the 
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sensor signal simultaneously with processing the first data set. The period corresponding to 
the first data set may overlap the period corresponding to the second data set. 

The control and measurement system may control the drive signal to maintain an 
amplitude of the sensor signal at a fixed setpoint, reduce the fixed setpoint when the drive 
5 signal exceeds a first threshold level, and increase the fixed setpoint when the drive signal is 
less than a second threshold level and the fixed setpoint is less than a maximum permitted 
value for the setpoint. The first threshold level may be 95% or less of a maximum permitted 
drive signal. 

The control and measurement system may perform an uncertainty analysis on the 
10 measurement. In this case, the control and measurement system may transmit the 
measurement and results of the uncertainty analysis to the control system. 

The control and measurement system may use digital processing to adjust a phase 
of the drive signal to compensate for a time delay associated with the sensor and components 
connected between the sensor and the driver. 
15 The details of one or more implementations are set forth in the accompanying 

drawings and the description below. Other features and advantages will be apparent from the 
description and drawings, and from the claims. 



DESCRIPTION OF DRAWINGS 

20 Fig. 1 is a block diagram of a digital mass flowmeter. 

Figs. 2 A and 2B are perspective and side views of mechanical components of a 
mass flowmeter. 

Figs. 3A-3C are schematic representations of three modes of motion of the 
flowmeter of Fig. 1. 

25 Fig. 4 is a block diagram of an analog control and measurement circuit. 

Fig. 5 is a block diagram of a digital mass flowmeter. 

Fig. 6 is a flow chart showing operation of the meter of Fig. 5. 

Figs. 7A and 7B are graphs of sensor data. 

Figs. 8A and 8B are graphs of sensor voltage relative to time. 
30 Fig. 9 is a flow chart of a curve fitting procedure. 

Fig. 10 is a flow chart of a procedure for generating phase differences. 
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Figs. 11A-11D, 12A-12D, and 13A-13D illustrate drive and sensor voltages at 
system startup. 

Fig. 14 is a flow chart of a procedure for measuring frequency, amplitude, and 
phase of sensor data using a synchronous modulation technique. 

Figs. 1 5 A and 1 5B are block diagrams of a mass flowmeter. 

Fig. 16 is a flow chart of a procedure implemented by the meter of Figs. 15A and 

15B. 

Fig. 17 illustrates log-amplitude control of a transfer function. 
Fig. 18 is a root locus diagram. 

Figs. 19A-19D are graphs of analog-to-digital converter performance relative to 
temperature. 

Figs. 20A-20C are graphs of phase measurements. 

Figs. 21A and 21B are graphs of phase measurements. 

Fig. 22 is a flow chart of a zero offset compensation procedure. 

Figs. 23A-23C, 24A, and 24B are graphs of phase measurements. 

Fig. 25 is a graph of sensor voltage. 

Fig. 26 is a flow chart of a procedure for compensating for dynamic effects. 
Figs. 27A-35E are graphs illustrating application of the procedure of Fig. 29. 
Figs. 36A-36L are graphs illustrating phase measurement. 
Fig. 37A is a graph of sensor voltages. 

Figs. 37B and 37C are graphs of phase and frequency measurements 
corresponding to the sensor voltages of Fig. 3 7 A. 

Figs. 37D and 37E are graphs of correction parameters for the phase and 
frequency measurements of Figs. 37B and 37C. 

Figs. 38A-38H are graphs of raw measurements. 

Figs. 39A-39H are graphs of corrected measurements. 

Figs. 40A-40H are graphs illustrating correction for aeration. 

Fig. 41 is a block diagram illustrating the effect of aeration in a conduit. 

Fig. 42 is a flow chart of a setpoint control procedure. 

Figs. 43A-43C are graphs illustrating application of the procedure of Fig. 41. 

Fig. 44 is a graph comparing the performance of digital and analog flowmeters. 
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Fig. 45 is a flow chart showing operation of a self-validating meter. 
Fig. 46 is a block diagram of a two-wire digital mass flowmeter. 

DETAILED DESCRIPTION 

5 Referring to Fig. 1, a digital mass flowmeter 100 includes a digital controller 105, 

one or more motion sensors 110, one or more drivers 1 15, a conduit 120 (also referred to as a 
flowtube), and a temperature sensor 125. The digital controller 105 may be implemented 
using one or more of, for example, a processor, a field-programmable gate array, an ASIC, 
other programmable logic or gate arrays, or programmable logic with a processor core. The 
10 digital controller generates a measurement of mass flow through the conduit 120 based at 
least on signals received from the motion sensors 110. The digital controller also controls the 
drivers 1 15 to induce motion in the conduit 120. This motion is sensed by the motion sensors 
110. 

Mass flow through the conduit 120 is related to the motion induced in the conduit 
15 in response to a driving force supplied by the drivers 115. In particular, mass flow is related 
to the phase and frequency of the motion, as well as to the temperature of the conduit. The 
digital mass flowmeter also may provide a measurement of the density of material flowing 
through the conduit. The density is related to the frequency of the motion and the 
temperature of the conduit. Many of the described techniques are applicable to a 
20 densitometer that provides a measure of density rather than a measure of mass flow. 

The temperature in the conduit, which is measured using the temperature sensor 
125, affects certain properties of the conduit, such as its stiffness and dimensions. The digital 
controller compensates for these temperature effects. The temperature of the digital 
controller 105 affects, for example, the operating frequency of the digital controller. In 
25 general, the effects of controller temperature are sufficiently small to be considered 

negligible. However, in some instances, the digital controller may measure the controller 
temperature using a solid state device and may compensate for effects of the controller 
temperature. 
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A. Mechanical Design 

In one implementation, as illustrated in Figs. 2 A and 2B, the conduit 120 is 
designed to be inserted in a pipeline (not shown) having a small section removed or reserved 
to make room for the conduit. The conduit 120 includes mounting flanges 12 for connection 
5 to the pipeline, and a central manifold block 16 supporting two parallel planar loops 18 and 
20 that are oriented perpendicularly to the pipeline. An electromagnetic driver 46 and a 
sensor 48 are attached between each end of loops 18 and 20. Each of the two drivers 46 
corresponds to a driver 115 of Fig. 1, while each of the two sensors 48 corresponds to a 
sensor 120 of Fig. 1. 

10 The drivers 46 on opposite ends of the loops are energized with current of equal 

magnitude but opposite sign (i.e., currents that are 180° out-of-phase) to cause straight 
sections 26 of the loops 18, 20 to rotate about their co-planar perpendicular bisector 56, 
which intersects the tube at point P (Fig. 2B). Repeatedly reversing (e.g., controlling 
sinusoidally) the energizing current supplied to the drivers causes each straight section 26 to 

15 undergo oscillatory motion that sweeps out a bow tie shape in the horizontal plane about line 
56-56, the axis of symmetry of the loop. The entire lateral excursion of the loops at the lower 
rounded turns 38 and 40 is small, on the order of 1/16 of an inch for a two foot long straight 
section 26 of a pipe having a one inch diameter. The frequency of oscillation is typically 
about 80 to 90 Hertz. 

20 

B. Conduit Motion 

The motion of the straight sections of loops 18 and 20 are shown in three modes 
in Figs. 3 A, 3B and 3C. In the drive mode shown in Fig. 3B, the loops are driven 180° out- 
of-phase about their respective points P so that the two loops rotate synchronously but in the 
25 opposite sense. Consequently, respective ends such as A and C periodically come together 
and go apart. 

The drive motion shown in Fig. 3B induces the Coriolis mode motion shown in 
Fig. 3 A, which is in opposite directions between the loops and moves the straight sections 26 
slightly toward (or away) from each other. The Coriolis effect is directly related to mvW, 
30 where m is the mass of material in a cross section of a loop, v is the velocity at which the 
mass is moving (the volumetric flow rate), W is the angular velocity of the loop (W = 
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Wosincot), and mv is the mass flow rate. The Coriolis effect is greatest when the two straight 
sections are driven sinusoidaliy and have a sinusoidally varying angular velocity. Under 
these conditions, the Coriolis effect is 90° out-of-phase with the drive signal. 

Fig. 3C shows an undesirable common mode motion that deflects the loops in the 
5 same direction. This type of motion might be produced by an axial vibration in the pipeline 
in the embodiment of Figs. 2 A and 2B because the loops are perpendicular to the pipeline. 

The type of oscillation shown in Fig. 3B is called the antisymmetrical mode, and 
the Coriolis mode of Fig. 3 A is called the symmetrical mode. The natural frequency of 
oscillation in the antisymmetrical mode is a function of the torsional resilience of the legs. 
10 Ordinarily the resonant frequency of the antisymmetrical mode for conduits of the shape 
shown in Figs. 2A and 2B is higher than the resonant frequency of the symmetrical mode. 
To reduce the noise sensitivity of the mass flow measurement, it is desirable to maximize the 
Coriolis force for a given mass flow rate. As noted above, the loops are driven at their 
resonant frequency, and the Coriolis force is directly related to the frequency at which the 
15 loops are oscillating (i.e., the angular velocity of the loops). Accordingly, the loops are 
driven in the antisymmetrical mode, which tends to have the higher resonant frequency. 

Other implementations may include different conduit designs. For example, a 
single loop or a straight tube section may be employed as the conduit. 

20 C Electronic Design 

The digital controller 105 determines the mass flow rate by processing signals 
produced by the sensors 48 (i.e., the motion sensors 110) located at opposite ends of the 
loops. The signal produced by each sensor includes a component corresponding to the 
relative velocity at which the loops are driven by a driver positioned next to the sensor and a 

25 component corresponding to the relative velocity of the loops due to Coriolis forces induced 
in the loops. The loops are driven in the antisymmetrical mode, so that the components of 
the sensor signals corresponding to drive velocity are equal in magnitude but opposite in 
sign. The resulting Coriolis force is in the symmetrical mode so that the components of the 
sensor signals corresponding to Coriolis velocity are equal in magnitude and sign. Thus, 

30 differencing the signals cancels out the Coriolis velocity components and results in a 

difference that is proportional to the drive velocity. Similarly, summing the signals cancels 
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out the drive velocity components and results in a sum that is proportional to the Coriolis 
velocity, which, in turn, is proportional to the Coriolis force. This sum then may be used to 
determine the mass flow rate. 

5 1. Analog Control System 

The digital mass flowmeter 100 provides considerable advantages over 
traditional, analog mass flowmeters. For use in later discussion, Fig. 4 illustrates an analog 
control system 400 of a traditional mass flowmeter. The sensors 48 each produce a voltage 
signal, with signal V A o being produced by sensor 48a and signal Vbo being produced by 
10 sensor 48b. Vao and Vbo correspond to the velocity of the loops relative to each other at the 
positions of the sensors. Prior to processing, signals V A o and V B o are amplified at respective 
input amplifiers 405 and 410 to produce signals Vai and Vbi. To correct for imbalances in 
the amplifiers and the sensors, input amplifier 410 has a variable gain that is controlled by a 
balance signal coming from a feedback loop that contains a synchronous demodulator 415 
1 5 and an integrator 420. 

At the output of amplifier 405, signal V A i is of the form: 

V A1 = V D sin cot + V c cos cot, 

and, at the output of amplifier 410, signal V B i is of the form: 

V B] = -V D sin cot + V c cos cot, 

20 where V D and Vc are, respectively, the drive voltage and the Coriolis voltage, and co is the 
drive mode angular frequency. 

Voltages V A i and Vbi are differenced by operational amplifier 425 to produce: 

v drv = v ai - v bi = 2V D sin cot, 
where V DRV corresponds to the drive motion and is used to power the drivers. In addition to 
25 powering the drivers, V DRV is supplied to a positive going zero crossing detector 430 that 
produces an output square wave Fdrv having a frequency corresponding to that of Vdrv = 
27cFdrv). F drv is used as the input to a digital phase locked loop circuit 435. Fdrv also is 
supplied to a processor 440. 

Voltages V A i and V B i are summed by operational amplifier 445 to produce: 
30 V C0R = V A , + V B , = 2V C cos cot, 
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where V CO r is related to the induced Coriolis motion. 

Vcor is supplied to a synchronous demodulator 450 that produces an output 
voltage V M that is directly proportional to mass by rejecting the components of Vcor that do 
not have the same frequency as, and are not in phase with, a gating signal Q. The phase 
5 locked loop circuit 435 produces Q, which is a quadrature reference signal that has the same 
frequency (co) as Vdrv and is 90° out of phase with V D rv (i.e., in phase with Vcor)- 
Accordingly, synchronous demodulator 450 rejects frequencies other than co so that V M 
corresponds to the amplitude of Vcor at co. This amplitude is directly proportional to the 
mass in the conduit. 

10 V M is supplied to a voltage-to-frequency converter 455 that produces a square 

wave signal Fm having a frequency that corresponds to the amplitude of Vm- The processor 
440 then divides F M by F D rv to produce a measurement of the mass flow rate. 

Digital phase locked loop circuit 435 also produces a reference signal I that is in 
phase with Vdrv and is used to gate the synchronous demodulator 415 in the feedback loop 

15 controlling amplifier 410. When the gains of the input amplifiers 405 and 410 multiplied by 
the drive components of the corresponding input signals are equal, the summing operation at 
operational amplifier 445 produces zero drive component (i.e., no signal in phase with Vdrv) 
in the signal Vcor- When the gains of the input amplifiers 405 and 410 are not equal, a drive 
component exists in V C or. This drive component is extracted by synchronous demodulator 

20 415 and integrated by integrator 420 to generate an error voltage that corrects the gain of 
input amplifier 410. When the gain is too high or too low, the synchronous demodulator 415 
produces an output voltage that causes the integrator to change the error voltage that modifies 
the gain. When the gain reaches the desired value, the output of the synchronous modulator 
goes to zero and the error voltage stops changing to maintain the gain at the desired value. 
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2. Digital Control System 

Fig. 5 provides a block diagram of an implementation 500 of the digital mass 
flowmeter 100 that includes the conduit 120, drivers 46, and sensors 48 of Figs. 2 A and 2B, 
along with a digital controller 505. Analog signals from the sensors 48 are converted to 
digital signals by analog-to-digital ("A/D") converters 510 and supplied to the controller 505. 
The A/D converters may be implemented as separate converters, or as separate channels of a 
single converter. 

Digital-to-analog ("D/A") converters 515 convert digital control signals from the 
controller 505 to analog signals for driving the drivers 46. The use of a separate drive signal 
for each driver has a number of advantages. For example, the system may easily switch 
between symmetrical and antisymmetrical drive modes for diagnostic purposes. In other 
implementations, the signals produced by converters 515 may be amplified by amplifiers 
prior to being supplied to the drivers 46. In still other implementations, a single D/A 
converter may be used to produce a drive signal applied to both drivers, with the drive signal 
being inverted prior to being provided to one of the drivers to drive the conduit 120 in the 

antisymmetrical mode. 

High precision resistors 520 and amplifiers 525 are used to measure the current 
supplied to each driver 46. A/D converters 530 convert the measured current to digital 
signals and supply the digital signals to controller 505. The controller 505 uses the measured 
currents in generating the driving signals. 

Temperature sensors 535 and pressure sensors 540 measure, respectively, the 
temperature and the pressure at the inlet 545 and the outlet 550 of the conduit. A/D 
converters 555 convert the measured values to digital signals and supply the digital signals to 
the controller 505. The controller 505 uses the measured values in a number of ways. For 
example, the difference between the pressure measurements may be used to determine a back 
pressure in the conduit. Since the stifmess of the conduit varies with the back pressure, the 
controller may account for conduit stiffness based on the determined back pressure. 

An additional temperature sensor 560 measures the temperature of the crystal 
oscillator 565 used by the A/D converters. An A/D converter 570 converts this temperature 
measurement to a digital signal for use by the controller 505. The input/output relationship 
of the A/D converters varies with the operating frequency of the converters, and the 
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operating frequency varies with the temperature of the crystal oscillator. Accordingly, the 
controller uses the temperature measurement to adjust the data provided by the A/D 
converters, or in system calibration. 

In the implementation of Fig. 5, the digital controller 505 processes the digitized 
5 sensor signals produced by the A/D converters 510 according to the procedure 600 illustrated 
in Fig. 6 to generate the mass flow measurement and the drive signal supplied to the drivers 
46. Initially, the controller collects data from the sensors (step 605). Using this data, the 
controller determines the frequency of the sensor signals (step 610), eliminates zero offset 
from the sensor signals (step 615), and determines the amplitude (step 620) and phase (step 

10 625) of the sensor signals. The controller uses these calculated values to generate the drive 
signal (step 630) and to generate the mass flow and other measurements (step 635). After 
generating the drive signals and measurements, the controller collects a new set of data and 
repeats the procedure. The steps of the procedure 600 may be performed serially or in 
parallel, and may be performed in varying order. 

15 Because of the relationships between frequency, zero offset, amplitude, and 

phase, an estimate of one may be used in calculating another. This leads to repeated 
calculations to improve accuracy. For example, an initial frequency determination used in 
determining the zero offset in the sensor signals may be revised using offset-eliminated 
sensor signals. In addition, where appropriate, values generated for a cycle may be used as 

20 starting estimates for a following cycle. 

Fig. 5 provides a general description of the hardware included in a digital 
flowmeter. A more detailed description of a specific hardware implementation is provided in 
the attached appendix, which is incorporated by reference. 

25 a. Data Collection 

For ease of discussion, the digitized signals from the two sensors will be referred 
to as signals SVi and SV2, with signal SVi coming from sensor 48a and signal SV2 coming 
from sensor 48b. Although new data is generated constantly, it is assumed that calculations 
are based upon data corresponding to one complete cycle of both sensors. With sufficient 

30 data buffering, this condition will be true so long as the average time to process data is less 
than the time taken to collect the data. Tasks to be carried out for a cycle include deciding 
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that the cycle has been completed, calculating the frequency of the cycle (or the frequencies 
of SVi and SV 2 ), calculating the amplitudes of SVi and SV 2 , and calculating the phase 
difference between SVi and SV2, In some implementations, these calculations are repeated 
for each cycle using the end point of the previous cycle as the start for the next. In other 

5 implementations, the cycles overlap by 180° or other amounts (e.g., 90°) so that a cycle is 
subsumed within the cycles that precede and follow it. 

Figs. 7A and 7B illustrate two vectors of sampled data from signals SVi and SV2, 
which are named, respectively, svl_jn and sv2_in. The first sampling point of each vector is 
known, and corresponds to a zero crossing of the sine wave represented by the vector. For 

10 svl_in, the first sampling point is the zero crossing from a negative value to a positive value, 
while for sv2_in the first sampling point is the zero crossing from a positive value to a 
negative value. 

An actual starting point for a cycle (i.e., the actual zero crossing) will rarely 
coincide exactly with a sampling point. For this reason, the initial sampling points 
15 (start_sample_SVl and start_sample_SV2) are the sampling points occurring just before the 
start of the cycle. To account for the difference between the first sampling point and the 
actual start of the cycle, the approach also uses the position (startoffsetSVl or 
start_offset_SV2) between the starting sample and the next sample at which the cycle 
actually begins. 

20 Since there is a phase offset between signals SV] and SV2, svl_in and sv2_in may 

start at different sampling points. If both the sample rate and the phase difference are high, 
there may be a difference of several samples between the start of svl_in and the start of 
sv2_in. This difference provides a crude estimate of the phase offset, and may be used as a 
check on the calculated phase offset, which is discussed below. For example, when sampling 

25 at 55 kHz, one sample corresponds to approximately 0.5 degrees of phase shift, and one cycle 
corresponds to about 800 sample points. 

When the controller employs functions such as the sum (A+B) and difference (A- 
B), with B weighted to have the same amplitude as A, additional variables (e.g., 
start_sample_sum and start_offset_sum) track the start of the period for each function. The 

30 sum and difference functions have a phase offset halfway between SVj and SV2. 
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In one implementation, the data structure employed to store the data from the 
sensors is a circular list for each sensor, with a capacity of at least twice the maximum 
number of samples in a cycle. With this data structure, processing may be carried out on data 
for a current cycle while interrupts or other techniques are used to add data for a following 
5 cycle to the lists. 

Processing is performed on data corresponding to a full cycle to avoid errors 
when employing approaches based on sine-waves. Accordingly, the first task in assembling 
data for a cycle is to determine where the cycle begins and ends. When nonoverlapping 
cycles are employed, the beginning of the cycle may be identified as the end of the previous 
10 cycle. When overlapping cycles are employed, and the cycles overlap by 180°, the beginning 
of the cycle may be identified as the midpoint of the previous cycle, or as the endpoint of the 
cycle preceding the previous cycle. 

The end of the cycle may be first estimated based on the parameters of the 
previous cycle and under the assumption that the parameters will not change by more than a 
15 predetermined amount from cycle to cycle. For example, five percent may be used as the 
maximum permitted change from the last cycle's value, which is reasonable since, at 
sampling rates of 55 kHz, repeated increases or decreases of five percent in amplitude or 
frequency over consecutive cycles would result in changes of close to 5,000 percent in one 
second. 

20 By designating five percent as the maximum permissible increase in amplitude 

and frequency, and allowing for a maximum phase change of 5° in consecutive cycles, a 
conservative estimate for the upper limit on the end of the cycle for signal SVi may be 
determined as: 



end_sample_SVl < start_sample_SVl + 



365 # sample_rate 



360 est_freq * 0.95 

25 where start_sampleJSVl is the first sample of svl_in, sample_rate is the sampling rate, and 
est_freq is the frequency from the previous cycle. The upper limit on the end of the cycle for 
signal SV 2 (end_sample_SV2) may be determined similarly. 

After the end of a cycle is identified, simple checks may be made as to whether 
the cycle is worth processing. A cycle may not be worth processing when, for example, the 
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conduit has stalled or the sensor waveforms are severely distorted. Processing only suitable 

cycles provides considerable reductions in computation. 

One way to determine cycle suitability is to examine certain points of a cycle to 

confirm expected behavior. As noted above, the amplitudes and frequency of the last cycle 
5 give useful starting estimates of the corresponding values for the current cycle. Using these 

values, the points corresponding to 30°, 150°, 210° and 330° of the cycle may be examined. 

If the amplitude and frequency were to match exactly the amplitude and frequency for the 

previous cycle, these points should have values corresponding to est_amp/2, est_amp/2, - 

est_amp/2, and -est_amp/2, respectively, where est_amp is the estimated amplitude of a 
10 signal (i.e., the amplitude from the previous cycle). Allowing for a five percent change in 

both amplitude and frequency, inequalities may be generated for each quarter cycle. For the 

30° point, the inequality is 



svl in 



^ 1 oui 30 * sample rate N 
start_sample_SVl + * 



> 0.475 * est_amp_S V 



360 est_freq*l.05j 

The inequalities for the other points have the same form, with the degree offset term (x/360) 

15 and the sign of the est_amp_SVl term having appropriate values. These inequalities can be 

used to check that the conduit has vibrated in a reasonable manner. 

Measurement processing takes place on the vectors svl Jn(start:end) and 

sv2_in(start:end) where: 

start = min (start_sample_SVl, start_sample_SV2), and 
end = max (end_sample_SVl, end_sample_SV2). 

20 

The difference between the start and end points for a signal is indicative of the frequency of 
the signal. 

b. Frequency Determination 
25 The frequency of a discretely-sampled pure sine wave may be calculated by 

detecting the transition between periods (i.e., by detecting positive or negative zero- 
crossings) and counting the number of samples in each period. Using this method, sampling, 
for example, an 82.2 Hz sine wave at 55 kHz will provide an estimate of frequency with a 
maximum error of 0.15 percent. Greater accuracy may be achieved by estimating the 
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fractional part of a sample at which the zero-crossing actually occurred using, for example, 
start_offsetjSVl and start_offset_SV2. Random noise and zero offset may reduce the 
accuracy of this approach. 

As illustrated in Figs. 8A and 8B, a more refined method of frequency 
determination uses quadratic interpolation of the square of the sine wave. With this method, 
the square of the sine wave is calculated, a quadratic function is fitted to match the minimum 
point of the squared sine wave, and the zeros of the quadratic function are used to determine 
the frequency. If 

sv t = A sin x t + 6 + as t 

where sv t is the sensor voltage at time t, A is the amplitude of oscillation, x t is the radian 
angle at time t (i.e., x t = 27ift), 5 is the zero offset, a t is a random variable with distribution 
N(0,1), and a is the variance of the noise, then the squared function is given by: 

sv 2 = A 2 sin 2 x, + 2 A (<5 + <T£ t ) sin * t + 28<7€ t + 5 2 + cr 2 £, 2 . 

When x t is close to 2tc, sin x t and sin 2 x t can be approximated as x ot = x r 2n and x 0t 2 , 
respectively. Accordingly, for values of x t close to 2%, a t can be approximated as: 

a 2 * A 2 x 2 t +2A (£ + os t )x 0t + 2Sae t +S 2 + <r 2 e 2 
» (A 2 * 2 t +2A<Sx 0t +£ 2 )+ cr£ t (2Ax 0t +2S + <7€ { ). 

This is a pure quadratic (with a non-zero minimum, assuming 5 7^0) plus noise, with the 
amplitude of the noise being dependent upon both a and 6. Linear interpolation also could 
be used. 

Error sources associated with this curve fitting technique are random noise, zero 
offset, and deviation from a true quadratic. Curve fitting is highly sensitive to the level of 
random noise. Zero offset in the sensor voltage increases the amplitude of noise in the sine- 
squared function, and illustrates the importance of zero offset elimination (discussed below). 
Moving away from the minimum, the square of even a pure sine wave is not entirely 
quadratic. The most significant extra term is of fourth order. By contrast, the most 
significant extra term for linear interpolation is of third order. 
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Degrees of freedom associated with this curve fitting technique are related to how 
many, and which, data points are used. The minimum is three, but more may be used (at 
greater computational expense) by using least-squares fitting. Such a fit is less susceptible to 
random noise. Fig. 8A illustrates that a quadratic approximation is good up to some 20° 

5 away from the minimum point. Using data points further away from the minimum will 
reduce the influence of random noise, but will increase the errors due to the non-quadratic 
terms (i.e., fourth order and higher) in the sine-squared function. 

Fig. 9 illustrates a procedure 900 for performing the curve fitting technique. As a 
first step, the controller initializes variables (step 905). These variables include end__point, 

10 the best estimate of the zero crossing point; epjnt, the integer value nearest to end_point; 
s[0,.i], the set of all sample points; z[k], the square of the sample point closest to endjoint; 
z[0..n-l], a set of squared sample points used to calculate end_point; n, the number of sample 
points used to calculate endjpoint (n=2k+l); step_length, the number of samples in s 
between consecutive values in z; and iteration_count, a count of the iterations that the 

1 5 controller has performed. 

The controller then generates a first estimate of end__point (step 910). The 
controller generates this estimate by calculating an estimated zero-crossing point based on 
the estimated frequency from the previous cycle and searching around the estimated crossing 
point (forwards and backwards) to find the nearest true crossing point (i.e., the occurrence of 

20 consecutive samples with different signs). The controller then sets end_point equal to the 
sample point having the smaller magnitude of the samples surrounding the true crossing 
point. 

Next, the controller sets n, the number of points for curve fitting (step 915). The 
controller sets n equal to 5 for a sample rate of 1 1 kHz, and to 21 for a sample rate of 44 kHz. 
25 The controller then sets iteration_count to 0 (step 920) and increments iteration_count (step 
925) to begin the iterative portion of the procedure. 

As a first step in the iterative portion of the procedure, the controller selects 
stepjength (step 930) based on the value of iteration_count. The controller sets step Jength 
equal to 6, 3, or 1 depending on whether iteration_count equals, respectively, 1, 2 or 3. 
30 Next, the controller determines ep _int as the integer portion of the sum of 

endjoint and 0.5 (step 935) and fills the z array (step 940). For example, when n equals 5, 
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z[0] = s[ep_int-2*stepjength] 2 , z[l] = s[ep_int-stepjength] 2 , z[2] = s[ep__int] 2 , 
z[3] = s[ep_int+stepjength] 2 , and z[4] = s[ep_int+2*stepjength] 2 . 

Next, the controller uses a filter, such as a Savitzky-Golay filter, to calculate 
smoothed values of z[k-l], z[k] and z[k+l] (step 945). Savitzky-Golay smoothing filters are 
5 discussed by Press et al. in Numerical Recipes in C pp. 650-655 (2nd ed., Cambridge 
University Press, 1995), which is incorporated by reference. The controller then fits a 

quadratic to z[k-l], z[k] and z[k+l] (step 950), and calculates the minimum value of the 

* * 

quadratic (z ) and the corresponding position (x ) (step 955). 

If x* is between the points corresponding to k-1 and k+1 (step 960), then the 
10 controller sets end_point equal to x* (step 965). Thereafter, if iteration_count is less than 3 

(step 970), the controller increments iteration_count (step 925) and repeats the iterative 

portion of the procedure. 

If x is not between the points corresponding to k-1 and k+1 (step 960), or if 

iteration_count equals 3 (step 970), the controller exits the iterative portion of the procedure. 
15 The controller then calculates the frequency based on the difference between end_point and 

the starting point for the cycle, which is known (step 975). 

In essence, the procedure 900 causes the controller to make three attempts to 

home in on end_point, using smaller step_lengths in each attempt. If the resulting minimum 

for any attempt falls outside of the points used to fit the curve (i.e., there has been 
20 extrapolation rather than interpolation), this indicates that either the previous or new estimate 

is poor, and that a reduction in step size is unwarranted. 

The procedure 900 may be applied to at least three different sine waves produced 

by the sensors. These include signals SVj and SV 2 and the weighted sum of the two. 

Moreover, assuming that zero offset is eliminated, the frequency estimates produced for these 
25 signals are independent. This is clearly true for signals SVi and SV2, as the errors on each 

are independent. It is also true, however, for the weighted sum, as long as the mass flow and 

the corresponding phase difference between signals SVi and SV 2 are large enough for the 

calculation of frequency to be based on different samples in each case. When this is true, the 

random errors in the frequency estimates also should be independent. 
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The three independent estimates of frequency can be combined to provide an 
improved estimate. This combined estimate is simply the mean of the three frequency 
estimates. 



5 c. Zero Offset Compensation 

An important error source in a Coriolis transmitter is zero offset in each of the 
sensor voltages. Zero offset is introduced into a sensor voltage signal by drift in the pre- 
amplification circuitry and the analog-to-digital converter. The zero offset effect may be 
worsened by slight differences in the pre-amplification gains for positive and negative 

10 voltages due to the use of differential circuitry. Each error source varies between 

transmitters, and will vary with transmitter temperature and more generally over time with 
component wear. 

An example of the zero offset compensation technique employed by the controller 
is discussed in detail below. In general, the controller uses the frequency estimate and an 
15 integration technique to determine the zero offset in each of the sensor signals. The 

controller then eliminates the zero offset from those signals. After eliminating zero offset 
from signals SVi and SV 2 , the controller may recalculate the frequency of those signals to 
provide an improved estimate of the frequency. 



20 d. Amplitude Determination 

The amplitude of oscillation has a variety of potential uses. These include 
regulating conduit oscillation via feedback, balancing contributions of sensor voltages when 
synthesizing driver waveforms, calculating sums and differences for phase measurement, and 
calculating an amplitude rate of change for measurement correction purposes. 

25 In one implementation, the controller uses the estimated amplitudes of signals 

SVi and SV 2 to calculate the sum and difference of signals SVi and SV 2 , and the product of 
the sum and difference. Prior to determining the sum and difference, the controller 
compensates one of the signals to account for differences between the gains of the two 
sensors. For example, the controller may compensate the data for signal S V 2 based on the 

30 ratio of the amplitude of signal SVj to the amplitude of signal SV 2 so that both signals have 
the same amplitude. 
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The controller may produce an additional estimate of the frequency based on the 
calculated sum. This estimate may be averaged with previous frequency estimates to 
produce a refined estimate of the frequency of the signals, or may replace the previous 
estimates. 

The controller may calculate the amplitude according to a Fourier-based 
technique to eliminate the effects of higher harmonics. A sensor voltage x(t) over a period T 
(as identified using zero crossing techniques) can be represented by an offset and a series of 
harmonic terms as: 

x(t) = V2 + aicos(oot) + a2cos(2cot) + a 3 cos(3cot) + ... + 
bisin(cot) + b 2 sin(2cot) + ... 

With this representation, a non-zero offset ao will result in non-zero cosine terms a n . Though 
the amplitude of interest is the amplitude of the fundamental component (i.e., the amplitude 
at frequency co), monitoring the amplitudes of higher harmonic components (i.e., at 
frequencies kco, where k is greater than 1) may be of value for diagnostic purposes. The 
values of a n and b n may be calculated as: 

2 CT 

and 

b = — f x(t)smnudt. 

J* JO 

The amplitude, A n , of each harmonic is given by: 

The integrals are calculated using Simpson's method with quadratic correction (described 
below). The chief computational expense of the method is calculating the pure sine and 
cosine functions. 

e. Phase Determination 
The controller may use a number of approaches to calculate the phase difference 
between signals SVi and SV2. For example, the controller may determine the phase offset of 
each harmonic, relative to the starting time at t=0, as: 



26 



Attorney Docket No.: 02052-079004 




The phase offset is interpreted in the context of a single waveform as being the difference 
between the start of the cycle (i.e., the zero-crossing point) and the point of zero phase for the 

5 component of SV(t) of frequency a>. Since the phase offset is an average over the entire 
waveform, it may be used as the phase offset from the midpoint of the cycle. Ideally, with no 
zero offset and constant amplitude of oscillation, the phase offset should be zero every cycle. 
The controller may determine the phase difference by comparing the phase offset of each 
sensor voltage over the same time period. 

10 The amplitude and phase may be generated using a Fourier method that eliminates 

the effects of higher harmonics. This method has the advantage that it does not assume that 
both ends of the conduits are oscillating at the same frequency. As a first step in the method, 
a frequency estimate is produced using the zero crossings to measure the time between the 
start and end of the cycle. If linear variation in frequency is assumed, this estimate equals the 

1 5 time-averaged frequency over the period. Using the estimated, and assumed time-invariant, 
frequency co of the cycle, the controller calculates: 



where SV(t) is the sensor voltage waveform (i.e., SVj(t) or SV 2 (t)). The controller then 
20 determines the estimates of the amplitude and phase: 




K 



In 



p" SV(t)cos(wt)dt, 



71 
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The controller then calculates a phase difference, assuming that the average phase 
and frequency of each sensor signal is representative of the entire waveform. Since these 
frequencies are different for SVi and SV2, the corresponding phases are scaled to the average 
frequency. In addition, the phases are shifted to the same starting point (i.e., the midpoint of 
5 the cycle on SVi). After scaling, they are subtracted to provide the phase difference: 

av_freq 



scaled_j)hase_SV 1 = phase_SVj 



frecLSV, 



scaled shift SV 2 = (midpoint,SV 2 - midpoint^ )h frec^Y, 
~ ~ 2 360 

scaled_phase_SV 2 =(phase_SV 2 + scale_shift_SV 2 ) av -^ ec * f 

freq_SV 2 



where h is the sample length and the midpoints are defined in terms of samples: 

(startpoint SV X + endpoint SV Y ) 
10 midpoint_SV x = - = — - = — — 

2 

In general, phase and amplitude are not calculated over the same time-frame for 
the two sensors. When the flow rate is zero, the two cycle mid-points are coincident. 
However, they diverge at high flow rates so that the calculations are based on sample sets 

15 that are not coincident in time. This leads to increased phase noise in conditions of changing 
mass flow. At full flow rate, a phase shift of 4° (out of 360°) means that only 99% of the 
samples in the SVi and SV2 data sets are coincident. Far greater phase shifts may be 
observed under aerated conditions, which may lead to even lower rates of overlap. 

Fig. 10 illustrates a modified approach 1000 that addresses this issue. First, the 

20 controller finds the frequencies (fj, f2) and the mid-points (mi, 1112) of the SVi and SV 2 data 
sets (di, d2) (step 1005). Assuming linear shift in frequency from the last cycle, the 
controller calculates the frequency of SV 2 at the midpoint of SVi (f2mi) and the frequency of 
SVi at the midpoint of SV 2 (fintf) (step 1010). 

The controller then calculates the starting and ending points of new data sets (d^n^ 

25 and d2mi) with mid-points m2 and mi respectively, and assuming frequencies of fim2 and f2mi 
(step 1015). These end points do not necessarily coincide with zero crossing points. 
However, this is not a requirement for Fourier-based calculations. 
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The controller then carries out the Fourier calculations of phase and amplitude on 
the sets di and d2mi, and the phase difference calculations outlined above (step 1020). Since 
the mid-points of di and d2mi are identical, scale-shift_SV 2 is always zero and can be ignored. 
The controller repeats these calculations for the data sets d 2 and d^ (step 1025). The 
5 controller then generates averages of the calculated amplitude and phase difference for use in 
measurement generation (step 1030). When there is sufficient separation between the mid 
points mi and m 2 , the controller also may use the two sets of results to provide local 
estimates of the rates of change of phase and amplitude. 

The controller also may use a difference-amplitude method that involves 

10 calculating a difference between SVi and SV 2 , squaring the calculated difference, and 

integrating the result. According to another approach, the controller synthesizes a sine wave, 
multiplies the sine wave by the difference between signals SVi and SV 2 , and integrates the 
result. The controller also may integrate the product of signals SVi and SV 2 , which is a sine 
wave having a frequency 2f (where f is the average frequency of signals SVi and SV 2 ), or 

15 may square the product and integrate the result. The controller also may synthesize a cosine 
wave comparable to the product sine wave and multiply the synthesized cosine wave by the 
product sine wave to produce a sine wave of frequency 4f that the controller then integrates. 
The controller also may use multiple ones of these approaches to produce separate phase 
measurements, and then may calculate a mean value of the separate measurements as the 

20 final phase measurement. 

The difference-amplitude method starts with: 



SV x (t) = A x sin 2rft + - and SV 2 (t) = A 2 sin 

v 2 y 



25 where cp is the phase difference between the sensors. Basic trigonometric identities may be 
used to define the sum (Sum) and difference (Diff) between the signals as: 

Sum = SV X (t) + SV 2 {i) = 2A X cos - • sin lift, and 

A 2 2 

Diff = SV X (0 - A SV 2 (t) = 2A X sin ^ * cos 2 itft. 

A 2 2 
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These functions have amplitudes of 2Aicos(<p2) and 2Aisin((p2), respectively. The controller 
calculates data sets for Sum and Diff from the data for SVj and SV2, and then uses one or 
more of the methods described above to calculate the amplitude of the signals represented by 
those data sets. The controller then uses the calculated amplitudes to calculate the phase 
5 difference, (p. 

As an alternative, the phase difference may be calculated using the function Prod, 

defined as: 

Prod = SumxDiff = 4A? cos — • sin — • cos 27ift • sin 27rft 

2 2 

= A x simp- sin 47rft , 

which is a function with amplitude A 2 sinq> and frequency 2f. Prod can be generated sample 
10 by sample, and cp may be calculated from the amplitude of the resulting sine wave. 

The calculation of phase is particularly dependent upon the accuracy of previous 
calculations (i.e., the calculation of the frequencies and amplitudes of SVi and SV2). The 
controller may use multiple methods to provide separate (if not entirely independent) 
estimates of the phase, which may be combined to give an improved estimate. 

15 

f . Drive Signal Generation 
The controller generates the drive signal by applying a gain to the difference 
between signals SVi and SV2. The controller may apply either a positive gain (resulting in 
positive feedback) or a negative gain (resulting in negative feedback). 
20 In general, the Q of the conduit is high enough that the conduit will resonate only 

at certain discrete frequencies. For example, the lowest resonant frequency for some 
conduits is between 65 Hz and 95 Hz, depending on the density of the process fluid, and 
irrespective of the drive frequency. As such, it is desirable to drive the conduit at the 
resonant frequency to minimize cycle-to-cycle energy loss. Feeding back the sensor voltage 
25 to the drivers permits the drive frequency to migrate to the resonant frequency. 

As an alternative to using feedback to generate the drive signal, pure sine waves 
having phases and frequencies determined as described above may be synthesized and sent to 
the drivers. This approach offers the advantage of eliminating undesirable high frequency 
components, such as harmonics of the resonant frequency. This approach also permits 



30 



Attorney Docket No.: 02052-079004 



compensation for time delays introduced by the analog-to-digital converters, processing, and 
digital-to-analog converters to ensure that the phase of the drive signal corresponds to the 
mid-point of the phases of the sensor signals. This compensation may be provided by 
determining the time delay of the system components and introducing a phase shift 
5 corresponding to the time delay. 

Another approach to driving the conduit is to use square wave pulses. This is 
another synthesis method, with fixed (positive and negative) direct current sources being 
switched on and off at timed intervals to provide the required energy. The switching is 
synchronized with the sensor voltage phase. Advantageously, this approach does not require 

10 digital-to-analog converters. 

In general, the amplitude of vibration of the conduit should rapidly achieve a 
desired value at startup, so as to quickly provide the measurement function, but should do so 
without significant overshoot, which may damage the meter. The desired rapid startup may 
be achieved by setting a very high gain so that the presence of random noise and the high Q 

15 of the conduit are sufficient to initiate motion of the conduit. In one implementation, high 
gain and positive feedback are used to initiate motion of the conduit. Once stable operation 
is attained, the system switches to a synthesis approach for generating the drive signals. 

Referring to Figs. 11A-13D, synthesis methods also may be used to initiate 
conduit motion when high gain is unable to do so. For example, if the DC voltage offset of 

20 the sensor voltages is significantly larger than random noise, the application of a high gain 
will not induce oscillatory motion. This condition is shown in Figs. 1 1 A- 1 ID, in which a 
high gain is applied at approximately 0.3 seconds. As shown in Figs. 1 1 A and 1 IB, 
application of the high gain causes one of the drive signals to assume a large positive value 
(Fig. 1 1 A) and the other to assume a large negative value (Fig. 1 IB). The magnitudes of the 

25 drive signals vary with noise in the sensor signals (Figs. 1 1C and 1 ID). However, the 

amplified noise is insufficient to vary the sign of the drive signals so as to induce oscillation. 

Figs. 12A-12D illustrate that imposition of a square wave over several cycles can 
reliably cause a rapid startup of oscillation. Oscillation of a conduit having a two inch 
diameter may be established in approximately two seconds. The establishment of conduit 

30 oscillation is indicated by the reduction in the amplitude of the drive signals, as shown in 



31 



Attorney Docket No.: 02052-079004 



Figs. 12A and 12B. Figs. 13A-13D illustrate that oscillation of a one inch conduit may be 
established in approximately half a second. 

A square wave also may be used during operation to correct conduit oscillation 
problems. For example, in some circumstances, flow meter conduits have been known to 
5 begin oscillating at harmonics of the resonant frequency of the conduit, such as frequencies 
on the order of 1.5 kHz. When such high frequency oscillations are detected, a square wave 
having a more desirable frequency may be used to return the conduit oscillation to the 
resonant frequency. 

10 g. Measurement Generation 

The controller digitally generates the mass flow measurement in a manner similar 
to the approach used by the analog controller. The controller also may generate other 
measurements, such as density. 

In one implementation, the controller calculates the mass flow based on the phase 
15 difference in degrees between the two sensor signals (phase_diff), the frequency of 
oscillation of the conduit (freq), and the process temperature (temp); 
T 2 = temp - T c , 

noneu_mf = tan (tc * phase_dif(7180), and 

massflow = 16 (MFi*T z 2 + MF 2 *T Z + MF 3 ) * noneumf/freq, 
20 where T c is a calibration temperature, MF r MF 3 are calibration constants calculated during a 
calibration procedure, and noneujnf is the mass flow in non-engineering units. 

The controller calculates the density based on the frequency of oscillation of the 
conduit and the process temperature: 

T z = temp - T c , 
25 c 2 = freq 2 , and 

density = (Di*T z 2 + D 2 *T Z + D 3 )/c 2 + D 4 *T Z 2 , 
where D1-D4 are calibration constants generated during a calibration procedure. 



D. Integration Techniques 
30 Many integration techniques are available, with different techniques requiring 

different levels of computational effort and providing different levels of accuracy. In the 
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described implementation, variants of Simpson's method are used. The basic technique may 
be expressed as: 



5 where tk is the time at sample k, yu is the corresponding function value, and h is the step 
length. This rule can be applied repeatedly to any data vector with an odd number of data 
points (i.e., three or more points), and is equivalent to fitting and integrating a cubic spline to 
the data points. If the number of data points happens to be even, then the so-called 3/8ths 
rule can be applied at one end of the interval: 



As stated earlier, each cycle begins and ends at some offset (e.g., 
start_offsetJSVl) from a sampling point. The accuracy of the integration techniques are 

15 improved considerably by taking these offsets into account. For example, in an integration of 
a half cycle sine wave, the areas corresponding to partial samples must be included in the 
calculations to avoid a consistent underestimate in the result. 

Two types of function are integrated in the described calculations: either sine or 
sine-squared functions. Both are easily approximated close to zero where the end points 

20 occur. At the end points, the sine wave is approximately linear and the sine-squared function 
is approximately quadratic. 

In view of these two types of functions, three different integration methods have 
been evaluated. These are Simpson's method with no end correction, Simpson's method with 
linear end correction, and Simpson's method with quadratic correction. 

25 The integration methods were tested by generating and sampling pure sine and 

sine-squared functions, without simulating any analog-to-digital truncation error. Integrals 
were calculated and the results were compared to the true amplitudes of the signals. The 
only source of error in these calculations was due to the integration techniques. The results 
obtained are illustrated in tables A and B. 




10 
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Table A - Integration of a sine function 



% err r (based on 1000 simutati ns) 


Av.Err rs (%) 


S.D. Error (%) 


Max.Error (%) 


Simpson Only 


-3.73e-3 


1 .33e-3 


6.17e-3 


Simpson + linear correction 


3.16e-8 


4.89e-8 


1 .56e-7 


Simpson + quadratic correction 


2.00e-4 


2.19e~2 


5.18e-1 



Table B - Integration of a sine-squared function 



■ 

% error (based on 1000 simulations 


Av.Errors (%) 


S.D. Error (%) 


Max.Error (%) 


Simpson Only 


-2.21 e-6 


1.10e-6 


4.39e-3 


Simpson + linear correction 


2.21 e-6 


6.93e-7 


2.52e-6 


Simpson + quadratic correction 


2.15e-11 


6.83e-11 


1.88e-10 



5 For sine functions, Simpson's method with linear correction was unbiased with 

the smallest standard deviation, while Simpson's method without correction was biased to a 
negative error and Simpson's method with quadratic correction had a relatively high standard 
deviation. For sine-squared functions, the errors were generally reduced, with the quadratic 
correction providing the best result. Based on these evaluations, linear correction is used 
10 when integrating sine functions and quadratic correction is used when integrating sine- 
squared functions. 

E. Synchronous Modulation Technique 

Fig. 14 illustrates an alternative procedure 1400 for processing the sensor signals. 

15 Procedure 1400 is based on synchronous modulation, such as is described by Denys et al., in 
"Measurement of Voltage Phase for the French Future Defence Plan Against Losses of 
Synchronism", IEEE Transactions on Power Delivery, 7(1), 62-69, 1992 and by Begovic et 
al. in "Frequency Tracking in Power Networks in the Presence of Harmonics", IEEE 
Transactions on Power Delivery, 8(2), 480-486, 1993, both of which are incorporated by 

20 reference. 

First, the controller generates an initial estimate of the nominal operating 
frequency of the system (step 1405). The controller then attempts to measure the deviation 
of the frequency of a signal x[k] (e.g., SVj) from this nominal frequency: 
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x[k] = A sin [(coo + Aco) kh + 3>] + e (£), 



10 



where A is the amplitude of the sine wave portion of the signal, coo is the nominal frequency 
(e.g., 88 Hz), Aco is the deviation from the nominal frequency, h is the sampling interval, <X> 
is the phase shift, and e(k) corresponds to the added noise and harmonics. 

To generate this measurement, the controller synthesizes two signals that oscillate 
at the nominal frequency (step 1410). The signals are phase shifted by 0 and 7i/2 and have 
amplitude of unity. The controller multiplies each of these signals by the original signal to 
produce signals yi and y2 (step 1415): 

^1 = x[k]cos( u 0 kh) = ysin[(2w 0 + Ao))kh + G>]+ y sin (Aco&/i +<X>),aw/ 

y 2 = x[k]sm(a) 0 kh) = ~ cos [(2^o + A«)tt + ®]+ ~ cos (Aw^A + <D), 

where the first terms of yi and y 2 are high frequency (e.g., 176 Hz) components and the 
second terms are low frequency (e.g., 0 Hz) components. The controller then eliminates the 
high frequency components using a low pass filter (step 1420): 

yj = Ysin(Aofc/i + ®)+ Ejfk], and 
y\ = yCos(AcoM +<D)+£ 2 [ k ]> 

where ei[k] and S2[k] represent the filtered noise from the original signals. The controller 
combines these signals to produce u[k] (step 1425): 

20 u[k] = (y([k] + jy' 2 [k])(y([k - 1] + jy 7 2 [k - 1]) 

= w 1 [k] + yw 2 [k] 

A 1 A 2 
= — cos (Aco h) + j — sin (Ao)h), 

4 4 

which carries the essential information about the frequency deviation. As shown, ui[k] 
25 represents the real component of u[k], while u 2 [k] represents the imaginary component. 

The controller uses the real and imaginary components of u[k] to calculate the 
frequency deviation, Af (step 1430): 



15 
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V =Iarctan^l. 
h «,[k] 

The controller then adds the frequency deviation to the nominal frequency (step 1435) to give 
the actual frequency: 

5 f = Af + f 0 . 

The controller also uses the real and imaginary components of u[k] to determine 
the amplitude of the original signal In particular, the controller determines the amplitude as 
(step 1440): 

10 A 1 =4<yjuf[k] + ul[k]. 



Next, the controller determines the phase difference between the two sensor 
signals (step 1445). Assuming that any noise (ei[k] and £2[k]) remaining after application of 
the low pass filter described below will be negligible, noise free versions of y\[k] and y2 ! [k] 
15 (yi*[k] and y2*[k]) may be expressed as: 

* A 
yj[k] = — sin(Acokh + O), and 

2 

* A 
y 2 [k] = — cos(Aco£/z - <£). 

2 

Multiplying these signals together gives: 



i2 



A 

v = y\y\ = — |sin(2<D) + sin(2Ao/:/2)| 
8 

20 Filtering this signal by a low pass filter having a cutoff frequency near 0 Hz removes the 
unwanted component and leaves: 

A 2 

v'=— sin(2<D), 

from which the phase difference can be calculated as: 

0 = — arcsm— . 
2 A 2 

25 This procedure relies on the accuracy with which the operating frequency is 

initially estimated, as the procedure measures only the deviation from this frequency. If a 
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good estimate is given, a very narrow filter can be used, which makes the procedure very 
accurate. For typical flowmeters, the operating frequencies are around 95 Hz (empty) and 82 
Hz (full). A first approximation of half range (88 Hz) is used, which allows a low-pass filter 
cut-off of 13 Hz. Care must be taken in selecting the cut-off frequency as a very small 
5 cut-off frequency can attenuate the amplitude of the sine wave. 

The accuracy of measurement also depends on the filtering characteristics 
employed. The attenuation of the filter in the dead-band determines the amount of harmonics 
rejection, while a smaller cutoff frequency improves the noise rejection. 

10 F. Meter with PI Control 

Figs. 15A and 15B illustrate a meter 1500 having a controller 1505 that uses 
another technique to generate the signals supplied to the drivers. Analog-to-digital 
converters 1510 digitize signals from the sensors 48 and provide the digitized signals to the 
controller 1505. The controller 1505 uses the digitized signals to calculate gains for each 

1 5 driver, with the gains being suitable for generating desired oscillations in the conduit. The 
gains may be either positive or negative. The controller 1505 then supplies the gains to 
multiplying digital-to-analog converters 1515. In other implementations, two or more 
multiplying digital-to-analog converters arranged in series may be used to implement a 
single, more sensitive multiplying digital-to-analog converter. 

20 The controller 1505 also generates drive signals using the digitized sensor signals. 

The controller 1505 provides these drive signals to digital-to-analog converters 1520 that 
convert the signals to analog signals that are supplied to the multiplying digital-to-analog 
converters 1515. 

The multiplying digital-to-analog converters 1515 multiply the analog signals by 
25 the gains from the controller 1505 to produce signals for driving the conduit. Amplifiers 
1525 then amplify these signals and supply them to the drivers 46. Similar results could be 
obtained by having the controller 1505 perform the multiplication performed by the 
multiplying digital-to-analog converter, at which point the multiplying digital-to-analog 
converter could be replaced by a standard digital-to-analog converter. 
30 Fig. 15B illustrates the control approach in more detail. Within the controller 

1505, the digitized sensor signals are provided to an amplitude detector 1550, which 
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determines a measure, a(t), of the amplitude of motion of the conduit using, for example, the 
technique described above. A summer 1555 then uses the amplitude a(t) and a desired 
amplitude ao to calculate an error e(t) as: 

e (t) = a 0 -a (t). 

5 The error e(t) is used by a proportional-integral ("PI") control block 1560 to generate a gain 
Ko(t). This gain is multiplied by the difference of the sensor signals to generate the drive 
signal. The PI control block permits high speed response to changing conditions. The 
amplitude detector 1550, summer 1555, and PI control block 1560 may be implemented as 
software processed by the controller 1505, or as separate circuitry. 

10 

1. Control Procedure 

The meter 1500 operates according to the procedure 1600 illustrated in Fig. 16. 
Initially, the controller receives digitized data from the sensors (step 1605). Thereafter, the 
procedure 1600 includes three parallel branches: a measurement branch 1610, a drive signal 

1 5 generation branch 1615, and a gain generation branch 1 620. 

In the measurement branch 1610, the digitized sensor data is used to generate 
measurements of amplitude, frequency, and phase, as described above (step 1625). These 
measurements then are used to calculate the mass flow rate (step 1630) and other process 
variables. In general, the controller 1505 implements the measurement branch 1610. 

20 In the drive signal generation branch 1615, the digitized signals from the two 

sensors are differenced to generate the signal (step 1635) that is multiplied by the gain to 
produce the drive signal. As described above, this differencing operation is performed by the 
controller 1505. In general, the differencing operation produces a weighted difference that 
accounts for amplitude differences between the sensor signals. 

25 In the gain generation branch 1620, the gain is calculated using the proportional- 

integral control block. As noted above, the amplitude, a(t), of motion of the conduit is 
determined (step 1640) and subtracted from the desired amplitude ao (step 1645) to calculate 
the error e(t). Though illustrated as a separate step, generation of the amplitude, a(t), may 
correspond to generation of the amplitude in the measurement generation step 1625. Finally, 

30 the PI control block uses the error e(t) to calculate the gain (step 1650). 
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The calculated gain is multiplied by the difference signal to generate the drive 
signal supplied to the drivers (step 1655). As described above, this multiplication operation 
is performed by the multiplying D/A converter or may be performed by the controller. 

5 2. PI Control Block 

The objective of the PI control block is to sustain in the conduit pure sinusoidal 
oscillations having an amplitude ao. The behavior of the conduit may be modeled as a simple 
mass-spring system that may be expressed as: 

x + 2C > (jy n x + ® 2 n x = 0 9 

10 where x is a function of time and the displacement of the mass from equilibrium, a) n is the 

natural frequency, and £ is a damping factor, which is assumed to be small (e.g., 0.001). 

The solution to this force equation as a function of an output y(t) and an input i(t) is 

analogous to an electrical network in which the transfer function between a supplied current, 

i(s), and a sensed output voltage, y(s), is: 

y(s) ks 

i(s) s 2 + 2(& n s + u 2 n 

To achieve the desired oscillation in the conduit, a positive-feedback loop having 
the gain Ko(t) is automatically adjusted by a 'slow' outer loop to give: 

x + (ity n - fcK Q (/)) x + 0) \x = 0. 

20 The system is assumed to have a "two-time-scales" property, which means that variations in 
Ko(t) are slow enough that solutions to the equation for x provided above can be obtained by 
assuming constant damping. 

A two-term PI control block that gives zero steady-state error may be expressed 

as: 

25 K,{t) = K p e{t) + K^e{t)dU 

where the error, e(t) (i.e., ao-a(t)), is the input to the PI control block, and K p and Ki are 
constants. In one implementation, with ao = 10, controller constants of K p =0.02 and 
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Kj=0.0005 provide a response in which oscillations build up quickly. However, this PI 
control block is nonlinear, which may result in design and operational difficulties. 

A linear model of the behavior of the oscillation amplitude may be derived by 
assuming that x(t) equals A&<°\ which results in: 

x = Ae jat +ja>e j(0t , and 

x = [A-co 2 A]e jo)t +2jcoAe jo)t . 

Substituting these expressions into the expression for oscillation of the loop, and separating 
into real and imaginary terms, gives: 

jco{2A + (2<;u n -kK 0 )A} = 0, and 

A + (2^ n -kK 0 )A + (c0n -cj 2 )a = 0. 

A(t) also may be expressed as: 

A kKn 
10 — = -Cw n + — -t. 

A 2 

A solution of this equation is: 



logA(t) = 



V 1 J 



t. 



15 Transforming variables by defining a(t) as being equal to log A(t), the equation for A(t) can 
be written as: 

da_ kK 0 (t) 

where Ko is now explicitly dependent on time. Taking Laplace transforms results in: 



a( S , = - kK » W/2 



which can be interpreted in terms of transfer-functions as in Fig. 17. This figure is of 
particular significance for the design of controllers, as it is linear for all Ko and a, with the 
only assumption being the two-time-scales property. The performance of the closed-loop is 
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10 



robust with respect to this assumption so that fast responses which are attainable in practice 
can be readily designed. 

From Fig. 17, the term £co n , is a "load disturbance" that needs to be eliminated by 
the controller (i.e., kKo/2 must equal £co n for a(t) to be constant). For zero steady-state error 
this implies that the outer-loop controller must have an integrator (or very large gain). As 
such, an appropriate PI controller, C(s), may be assumed to be Kp(l+l/sTj), where Ti is a 
constant. The proportional term is required for stability. The term £© n , however, does not 
affect stability or controller design, which is based instead on the open-loop transfer function: 

r/.w* a & *M + j3r '> **,/2(,-H/r,) 
= <^ = 2.T, = ? • 



The root locus for varying K p is shown in Fig. 18. For small K p , there are slow 
underdamped roots. As K p increases, the roots become real at the point P for which the 
controller gain is K p =8/(kTj). Note in particular that the theory does not place any restriction 
upon the choice of Ti. Hence the response can, in principle, be made critically damped and 
15 as fast as desired by appropriate choices of Kp and Ti. 

Although the poles are purely real at point P, this does not mean there is no 
overshoot in the closed-loop step response. This is most easily seen by inspecting the 
transfer function between the desired value, ao, and the error e: 

e(s) £^ s 2 

a 0 (s)~ s 2 +0.5^(5 + 1/^) ~ p 2 (sY 

20 where p 2 is a second-order polynomial. With a step input, ao(s)=a/s, the response can be 
written as ap'(t), where p(t) is the inverse transform of l/p2(s) and equals 
aiexp(-^it)+a 2 exp(-A,2t). The signal p(t) increases and then decays to zero so that e(t), which 
is proportional to p 1 , must change sign, implying overshoot in a(t). The set-point ao may be 
prefiltered to give a pseudo set-point a<,*: 

25 a 0 (s) = —^—a Q (s), 

1 + sT;. 

where Tj is the known controller parameter. With this prefilter, real controller poles should 
provide overshoot-free step responses. This feature is useful as there may be physical 
constraints on overshoot (e.g., mechanical interference or overstressing of components). 
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The root locus of Fig. 1 8 assumes that the only dynamics are from the inner-loop's 
gain/log-amplitude transfer function (Fig, 16) and the outer-loop's PI controller C(s) (i.e., that 
the log-amplitude a = log A is measured instantaneously). However, A is the amplitude of an 
oscillation which might be growing or decaying and hence cannot in general be measured 

5 without taking into account the underlying sinusoid. There are several possible methods for 
measuring A, in addition to those discussed above. Some are more suitable for use in 
quasi-steady conditions. For example, a phase-locked loop in which a sinusoidal signal 
s(t) = sin(cD n t+<J>o) locks onto the measured waveform y(t) = A(t)sin(a> n t+<l>i) may be 
employed. Thus, a measure of the amplitude a = log A is given by dividing these signals 

10 (with appropriate safeguards and filters). This method is perhaps satisfactory near the 
steady-state but not for start-up conditions before there is a lock. 

Another approach uses a peak-follower that includes a zero-crossing detector 
together with a peak- following algorithm implemented in the controller. Zero-crossing 
methods, however, can be susceptible to noise. In addition, results from a peak-follower are 

15 available only every half-cycle and thereby dictate the sample interval for controller updates. 

Finally, an AM detector may be employed. Given a sine wave y(t) = Asinco n t, an 
estimate of A may be obtained from Ai0.57cF{abs(y)}, where F{} is a suitable low-pass filter 
with unity DC gain. The AM detector is the simplest approach. Moreover, it does not 
presume that there are oscillations of any particular frequency, and hence is usable during 

20 startup conditions. It suffers from a disadvantage that there is a leakage of harmonics into 
the inner loop which will affect the spectrum of the resultant oscillations. In addition, the 
filter adds extra dynamics into the outer loop such that compromises need to be made 
between speed of response and spectral purity. In particular, an effect of the filter is to 
constrain the choice of the best Tj. 

25 The Fourier series for abs(y) is known to be: 



Aabs(sinco n t) = ^ 



2 2 2 

1 + — cos 2co n t cos 4co n t + — cos 6co n t + .. 

7T l 3 15 "35 n 

As such, the output has to be scaled by nil to give the correct DC output A, and the (even) 
harmonic terms akCOs2ka> n t have to be filtered out. As all the filter needs to do is to pass the 
DC component through and reduce all other frequencies, a "brick-wall" filter with a cut-off 
30 below 2o) n is sufficient. However, the dynamics of the filter will affect the behavior of the 
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closed-loop. A common choice of filter is in the Butterworth form. For example, the 
third-order low-pass filter with a design break-point frequency Ob is: 

F(s) = \ — r — — . 

l + 2s/u b + 2s /o)b + s /coj. 

At the design frequency the response is 3dB down; at 2cob it is -18dB (0.12), and at 4©b it is 
5 -36dB (0.01 5) down. Higher-order Butterworth filters have a steeper roll-off, but most of 
their poles are complex and may affect negatively the control-loop's root locus. 

G. Zero Offset Compensation 

As noted above, zero offset may be introduced into a sensor voltage signal by 

10 drift in the pre-amplification circuitry and by the analog-to-digital converter. Slight 

differences in the pre-amplification gains for positive and negative voltages due to the use of 
differential circuitry may worsen the zero offset effect. The errors vary between transmitters, 
and with transmitter temperature and component wear. 

Audio quality (i.e., relatively low cost) analog-to-digital converters may be 

15 employed for economic reasons. These devices are not designed with DC offset and 
amplitude stability as high priorities. Figs. 19A-19D show how offset and positive and 
negative gains vary with chip operating temperature for one such converter (the AD 1879 
converter). The repeatability of the illustrated trends is poor, and even allowing for 
temperature compensation based on the trends, residual zero offset and positive/negative gain 

20 mismatch remain. 

If phase is calculated using the time difference between zero crossing points on 
the two sensor voltages, DC offset may lead to phase errors. This effect is illustrated by Figs 
20A-20C. Each graph shows the calculated phase offset as measured by the digital 
transmitter when the true phase offset is zero (i.e., at zero flow). 

25 Fig. 20A shows phase calculated based on whole cycles starting with positive 

zero-crossings. The mean value is 0.00627 degrees. 

Fig. 20B shows phase calculated starting with negative zero-crossings. The mean 
value is 0.0109 degrees. 
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Fig. 20C shows phase calculated every half-cycle. Fig. 20C interleaves the data 
from Figs. 20A and 20B. The average phase (-0.00234) is closer to zero than in Figs. 20A 
and 20B, but the standard deviation of the signal is about six times higher. 

More sophisticated phase measurement techniques, such as those based on 
5 Fourier methods, are immune to DC offset. However, it is desirable to eliminate zero offset 
even when those techniques are used, since data is processed in whole-cycle packets 
delineated by zero crossing points. This allows simpler analysis of the effects of, for 
example, amplitude modulation on apparent phase and frequency. In addition, gain 
mismatch between positive and negative voltages will introduce errors into any measurement 
10 technique. 

The zero-crossing technique of phase detection may be used to demonstrate the 
impact of zero offset and gain mismatch error, and their consequent removal. Figs. 21 A and 
2 IB illustrate the long term drift in phase with zero flow. Each point represents an average 
over one minute of live data. Fig. 21 A shows the average phase, and Fig. 2 IB shows the 
15 standard deviation in phase. Over several hours, the drift is significant. Thus, even if the 
meter were zeroed every day, which in many applications would be considered an excessive 
maintenance requirement, there would still be considerable phase drift. 

1. Compensation Technique 
20 A technique for dealing with voltage offset and gain mismatch uses the 

computational capabilities of the digital transmitter and does not require a zero flow 
condition. The technique uses a set of calculations each cycle which, when averaged over a 
reasonable period (e.g., 10,000 cycles), and excluding regions of major change (e.g., set point 
change, onset of aeration), converge on the desired zero offset and gain mismatch 
25 compensations. 

Assuming the presence of up to three higher harmonics, the desired waveform for 
a sensor voltage SV(t) is of the form: 

SV(t) = A] sin(o?t) + A 2 sin(2o>t) + A 3 sin(3wt) + A 4 sin(4ut) 

where Ai designates the amplitude of the fundamental frequency component and A2-A4 
30 designate the amplitudes of the three harmonic components. However, in practice, the actual 
waveform is adulterated with zero offset Z 0 (which has a value close to zero) and mismatch 
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between the negative and positive gains G n and G p . Without any loss of generality, it can be 
assumed that G p equals one and that G n is given by: 

G n = 1 + e G , 
where 6g represents the gain mismatch. 

5 The technique assumes that the amplitudes Ai and the frequency o are constant. 

This is justified because estimates of Z 0 and Eq are based on averages taken over many 
cycles (e.g., 10,000 interleaved cycles occurring in about 1 minute of operation). When 
implementing the technique, the controller tests for the presence of significant changes in 
frequency and amplitude to ensure the validity of the analysis. The presence of the higher 

10 harmonics leads to the use of Fourier techniques for extracting phase and amplitude 
information for specific harmonics. This entails integrating SV(t) and multiplying by a 
modulating sine or cosine function. 

The zero offset impacts the integral limits, as well as the functional form. 
Because there is a zero offset, the starting point for calculation of amplitude and phase will 

15 not be at the zero phase point of the periodic waveform SV(t). For zero offset Z 0 , the 
corresponding phase offset is, approximately, 



20 



<Pz 0 = ~ sin 



r 7 \ 



For small phase, Ai with corresponding time delay 0 w 



The integrals are scaled so that the limiting value (i.e., as Z 0 and coo approach 
zero) equals the amplitude of the relevant harmonic. The first two integrals of interest are: 

lift = — ° f (SV(t) + Z 0 ).sin[«(t-t z Q )]dt, and 

71 



2 ~ +tz o 

2© 

MNs 



a z 0 

(1 + s G ) J (SV(t) + Z 0 ) • sin[w(t - t Zo )]dt. 
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10 



15 



These integrals represent what in practice is calculated during a normal Fourier analysis of 
the sensor voltage data. The subscript 1 indicates the first harmonic, N and P indicate, 
respectively, the negative or positive half cycle, and s and c indicate, respectively, whether a 
sine or a cosine modulating function has been used. 

Strictly speaking, the mid-zero crossing point, and hence the corresponding 
integral limits, should be given by n/co - tz 0> rather than 7i/co + tzo. However, the use of the 
exact mid-point rather than the exact zero crossing point leads to an easier analysis, and 
better numerical behavior (due principally to errors in the location of the zero crossing point). 
The only error introduced by using the exact mid-point is that a small section of each of the 
above integrals is multiplied by the wrong gain (1 instead of 1 + £q and vice versa). 
However, these errors are of order Z 0 6g and are considered negligible. 

Using computer algebra and assuming small Z 0 and Eg, first order estimates for 
the integrals may be derived as: 
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7T 
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hNs est - (1 + s g) 
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A x — Z 0 

7T 
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Useful related functions including the sum, difference, and ratio of the integrals 
and their estimates may be determined. The sum of the integrals may be expressed as: 

Sum ls = (/ 1A +I m ), 



while the sum of the estimates equals: 



Sum ls esf=4( 2 + s G) — Z 0 e G 

— 7T 



\ 2 ^ 4 4 

1 + — - + 

3 /\ 15 Ay 



20 Similarly, the difference of the integrals may be expressed as: 

Diff,, = 



while the difference of the estimates is: 



4 / x 
D *ff\$ est = A i e G +-Z 0 (2 + £ G ) 

TT 



3 A, 15 A, 



Finally, the ratio of the integrals is: 
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Ratio Xs = 



I 



IPs 
\Ns 



while the ratio of the estimates is: 



Ratio \ s est - 



-™ 1 + e. 



1 + Z 



o 



8 15^ + KM 2 + 4 ^4 



15 



Corresponding cosine integrals are defined as: 



10 



hpc - 



with estimates: 



7T J 'Z 0 



7T 



2c0 f l + ' Zo (SV(t) + Z 0 )cos[u(t-t 20 )]dt, and 

j jr + ' z °(5no+^o) co 4 w (^^ 0 )] dt ' 
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15 



and sums: 



Sum ic =I lPc +I lNc ,and 
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Sum lc est = e G 
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Second harmonic integrals are: 



I 2ft = — J f' z ° (SV(t) + Z 0 ) sin [2u(t - t Z0 )] dt, and 



7T * 1 zq 
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with estimates: 



J 2H est M 2 T , r ^0 



15* 



-5 + 9^ 
A 



i j 



,and 



est =0 + s c ) 
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and sums: 



Sum 2s =l 2 p s +hN S , and 
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Suni 2p, «, = A 2 (2 + e G )-— e G Z 0 
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The integrals can be calculated numerically every cycle. As discussed below, the 
equations estimating the values of the integrals in terms of various amplitudes and the zero 
offset and gain values are rearranged to give estimates of the zero offset and gain terms based 
on the calculated integrals. 



10 2. Example 

The accuracy of the estimation equations may be illustrated with an example. For 
each basic integral, three values are provided: the "true" value of the integral (calculated 
within Mathcad using Romberg integration), the value using the estimation equation, and the 
value calculated by the digital transmitter operating in simulation mode, using Simpson's 
1 5 method with end correction. 

Thus, for example, the value for Ii Ps calculated according to: 

A = — J F' Z ° ( SV « + z o)sin[o>(t - t ZQ )]dt 

is 0.101353, while the estimated value (Iip s est) calculated as: 

, 2 A 2 4 A 4 " 
1 + - + - 

3 A } 15 A x _ 

20 is 0. 101358. The value calculated using the digital transmitter in simulation mode is 
0.101340. These calculations use the parameter values illustrated in Table C. 
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Table C 



Parameter 


Value 


Comment 


(0 


1607T 


This corresponds to frequency = 80Hz, a typical value. 


A, 


0.1 


This more typically is 0.3, but it could be smaller with aeration. 


A 2 


0.01 


This more typically is 0.005, but it could be larger with aeration. 


A 3 and A 4 


0.0 


The digital Coriolis simulation mode only offers two harmonics, so these 
higher harmonics are ignored. However, they are small (<0.002). 


Z 0 


0.001 


Experience suggests that this is a large value for zero offset. 




0.001 


Experience suggests this is a large value for gain mismatch. 



The exact, estimate and simulation results from using these parameter values are illustrated in 
Table D. 
5 Table D 



Integral 


'Exact' Value 


Estimate 


Digital Coriolis Simulation 


I IPs 


0.101353 


0.101358 


0.101340 


IlNs 


0.098735 


0.098740 


0.098751 


IlPc 


0.007487 


0.007488 


0.007500 


IlNc 


-0.009496 


-0.009498 


-0.00953 1 


hps 


0.009149 


0.009151 


0.009118 


l2Ns 


0.010857 


0.010859 


0.010885 



Thus, at least for the particular values selected, the estimates given by the first 
order equations are extremely accurate. As Z 0 and e G approach zero, the errors in both the 
estimate and the simulation approach zero. 

10 

3. Implementation 

The first order estimates for the integrals define a series of non-linear equations in 
terms of the amplitudes of the harmonics, the zero offset, and the gain mismatch. As the 
equations are non-linear, an exact solution is not readily available. However, an 
15 approximation followed by corrective iterations provides reasonable convergence with 
limited computational overhead. 
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Conduit-specific ratios may be assumed for A1-A4. As such, no attempt is made 
to calculate all of the amplitudes A1-A4. Instead, only Ai and A2 are estimated using the 
integral equations defined above. Based on experience of the relative amplitudes, A3 may be 
approximated as A2/2, and A4 may be approximated as A2/IO. 
5 The zero offset compensation technique may be implemented according to the 

procedure 2200 illustrated in Fig. 22. During each cycle, the controller calculates the 
integrals Ii Ps , Ii Ns > Iipo Iino l2Ps> l2Ns and related functions sumi s , ration, sumi c and sum 2s 
(step 2205). This requires minimal additional calculation beyond the conventional Fourier 
calculations used to determine frequency, amplitude and phase. 
10 Every 10,000 cycles, the controller checks on the slope of the sensor voltage 

amplitude At, using a conventional rate-of-change estimation technique (step 2210). If the 
amplitude is constant (step 2215), then the controller proceeds with calculations for zero 
offset and gain mismatch. This check may be extended to test for frequency stability. 

To perform the calculations, the controller generates average values for the 
15 functions (e.g., sumi s ) over the last 10,000 cycles. The controller then makes a first 
estimation of zero offset and gain mismatch (step 2225): 

Zo = -Sumi c /2, and 
£g = 1/Ratiois - 1 

Using these values, the controller calculates an inverse gain factor (k) and 
20 amplitude factor (amp_factor) (step 2230): 

k= 1.0/(1.0 + 0.5 * e G ),and 
amp_factor = 1 + 50/75 * SunWSumis 
The controller uses the inverse gain factor and amplitude factor to make a first estimation of 
the amplitudes (step 2235): 
25 Ai = k * [Sumi s /2 + 2/n * Z 0 * £g * amp_factor], and 

A 2 = k * [Sum 2s /2 - 4/(3 * n) * Z 0 * e G 
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The controller then improves the estimate by the following calculations, iterating 
as required (step 2240): 

xi = Z 0 , 
x 2 = e Gj 

5 e G = [1 + 8/7t * Xj/At * amp_factor]/Ratioi s - 1.0, 

Z 0 = - Sumi c /2 + x 2 * (x, + 2.773/ti* A 2 )/2, 
Ai = k * [Sumi s /2 + 2/n * xi * x 2 * amp_factor], 
A 2 = k * [Sum 2s /2 - 4/(15 * n) * xi * x 2 * (5 - 4.5 * A 2 )]. 
The controller uses standard techniques to test for convergence of the values of Z 0 and e G . 
10 In practice the corrections are small after the first iteration, and experience suggests that three 
iterations are adequate. 

Finally, the controller adjusts the raw data to eliminate Z 0 and s G (step 2245). 
The controller then repeats the procedure. Once zero offset and gain mismatch have been 
eliminated from the raw data, the functions (i.e., sumi s ) used in generating subsequent values 
15 for Z 0 and e G are based on corrected data. Accordingly, these subsequent values for Z 0 and 
8 G reflect residual zero offset and gain mismatch, and are summed with previously generated 
values to produce the actual zero offset and gain mismatch. In one approach to adjusting the 
raw data, the controller generates adjustment parameters (e.g., Sl_off and S2_off) that are 
used in converting the analog signals from the sensors to digital data. 
20 Figs. 23 A-23C, 24A and 24B show results obtained using the procedure 2200. 

The short-term behavior is illustrated in Figs. 23A-23C. This shows consecutive phase 
estimates obtained five minutes after startup to allow time for the procedure to begin 
affecting the output. Phase is shown based on positive zero-crossings, negative zero- 
crossings, and both. 

25 The difference between the positive and negative mean values has been reduced 

by a factor of 20, with a corresponding reduction in mean zero offset in the interleaved data 
set. The corresponding standard deviation has been reduced by a factor of approximately 6. 

Longer term behavior is shown in Figs. 24A and 24B. The initial large zero offset 
is rapidly corrected, and then the phase offset is kept close to zero over many hours. The 

30 average phase offset, excluding the first few values, is 6.14e" 6 , which strongly suggests that 
the procedure is successful in compensating for changes in voltage offset and gain imbalance. 
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Typical values for Z G and €g for the digital Coriolis meter are Zn = -7.923c- 4 and 
e G = -1.754e 5 for signal SVi, and Z G = -8.038c- 4 and e G = +6.93Q 4 for signal SV 2 . 

H. Dynamic Analysis 

5 In general, conventional measurement calculations for Coriolis meters assume 

that the frequency and amplitude of oscillation on each side of the conduit are constant, and 
that the frequency on each side of the conduit is identical and equal to the so-called resonant 
frequency. Phases generally are not measured separately for each side of the conduit, and the 
phase difference between the two sides is assumed to be constant for the duration of the 

10 measurement process. Precise measurements of frequency, phase and amplitude every half- 
cycle using the digital meter demonstrate that these assumptions are only valid when 
parameter values are averaged over a time period on the order of seconds. Viewed at 100 Hz 
or higher frequencies, these parameters exhibit considerable variation. For example, during 
normal operation, the frequency and amplitude values of SVi may exhibit strong negative 

15 correlation with the corresponding SV 2 values. Accordingly, conventional measurement 
algorithms are subject to noise attributable to these dynamic variations. The noise becomes 
more significant as the measurement calculation rate increases. Other noise terms may be 
introduced by physical factors, such as flowtube dynamics, dynamic non-linearities (e.g. 
flowtube stiffness varying with amplitude), or the dynamic consequences of the sensor 

20 voltages providing velocity data rather than absolute position data. 

The described techniques exploit the high precision of the digital meter to monitor 
and compensate for dynamic conduit behavior to reduce noise so as to provide more precise 
measurements of process variables such as mass flow and density. This is achieved by 
monitoring and compensating for such effects as the rates of change of frequency, phase and 

25 amplitude, flowtube dynamics, and dynamic physical non-idealities. A phase difference 
calculation which does not assume the same frequency on each side has already been 
described above. Other compensation techniques are described below. 

Monitoring and compensation for dynamic effects may take place at the 
individual sensor level to provide corrected estimates of phase, frequency, amplitude or other 

30 parameters. Further compensation may also take place at the conduit level, where data from 
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both sensors are combined, for example in the calculation of phase difference and average 
frequency. These two levels may be used together to provide comprehensive compensation. 

Thus, instantaneous mass flow and density measurements by the flowmeter may 
be improved by modeling and accounting for dynamic effects of flowmeter operation. In 
5 general, 80% or more of phase noise in a Coriolis flowmeter may be attributed to flowtube 
dynamics (sometimes referred to as "ringing"), rather than to process conditions being 
measured. The application of a dynamic model can reduce phase noise by a factor of 4 to 10, 
leading to significantly improved flow measurement performance. A single model is 
effective for all flow rates and amplitudes of oscillation. Generally, computational 

10 requirements are negligible. 

The dynamic analysis may be performed on each of the sensor signals in isolation 
from the other. This avoids, or at least delays, modeling the dynamic interaction between the 
two sides of the conduit, which is likely to be far more complex than the dynamics at each 
sensor. Also, analyzing the individual sensor signals is more likely to be successful in 

15 circumstances such as batch startup and aeration where the two sides of the conduit are 
subject to different forces from the process fluid. 

In general, the dynamic analysis considers the impact of time- varying amplitude, 
frequency and phase on the calculated values for these parameters. While the frequency and 
amplitude are easily defined for the individual sensor voltages, phase is conventionally 

20 defined in terms of the difference between the sensor voltages. However, when a Fourier 
analysis is used, phase for the individual sensor may be defined in terms of the difference 
between the midpoint of the cycle and the average 180° phase point. 

Three types of dynamic effects are measurement error and the so-called 
"feedback" and "velocity" effects. Measurement error results because the algorithms for 

25 calculating amplitude and phase assume that frequency, amplitude, and phase are constant 
over the time interval of interest. Performance of the measurement algorithms may be 
improved by correcting for variations in these parameters. 

The feedback effect results from supplying energy to the conduit to make up for 
energy loss from the conduit so as to maintain a constant amplitude of oscillation. The need 

30 to add energy to the conduit is only recognized after the amplitude of oscillation begins to 
deviate from a desired setpoint. As a result, the damping term in the equation of motion for 
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the oscillating conduit is not zero, and, instead, constantly dithers around zero. Although the 
natural frequency of the conduit does not change, it is obscured by shifts in the zero- 
crossings (i.e., phase variations) associated with these small changes in amplitude. 

The velocity effect results because the sensor voltages observe conduit velocity, 
5 but are analyzed as being representative of conduit position. A consequence of this is that 
the rate of change of amplitude has an impact on the apparent frequency and phase, even if 
the true values of these parameters are constant. 



1. Sensor-Level Compensation for Amplitude Modulation 

10 One approach to correcting for dynamic effects monitors the amplitudes of the 

sensor signals and makes adjustments based on variations in the amplitudes. For purposes of 
analyzing dynamic effects, it is assumed that estimates of phase, frequency and amplitude 
may be determined for each sensor voltage during each cycle. As shown in Fig. 25, 
calculations are based on complete but overlapping cycles. Each cycle starts at a zero 

15 crossing point, halfway through the previous cycle. Positive cycles begin with positive 
voltages immediately after the initial zero-crossing, while negative cycles begin with 
negative voltages. Thus cycle n is positive, while cycles n-1 and n+ 1 are negative. It is 
assumed that zero offset correction has been performed so that zero offset is negligible. It 
also is assumed that higher harmonics may be present. 

20 Linear variation in amplitude, frequency, and phase are assumed. Under this 

assumption, the average value of each parameter during a cycle equals the instantaneous 
value of the parameter at the mid-point of the cycle. Since the cycles overlap by 180 
degrees, the average value for a cycle equals the starting value for the next cycle. 

For example, cycle n is from time 0 to 2tc/co. The average values of amplitude, 

25 frequency and phase equal the instantaneous values at the mid-point, 7i/co, which is also the 
starting point for cycle n+1, which is from time ti/co to 3n/a>. Of course, these timings are 
approximate, since co also varies with time. 



a. Dynamic Effect Compensation Procedure 
30 The controller accounts for dynamic effects according to the procedure 2600 

illustrated in Fig. 26. First, the controller produces a frequency estimate (step 2605) by using 
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the zero crossings to measure the time between the start and end of the cycle, as described 
above. Assuming that frequency varies linearly, this estimate equals the time-averaged 
frequency over the period. 

The controller then uses the estimated frequency to generate a first estimate of 
5 amplitude and phase using the Fourier method described above (step 2610). As noted above, 
this method eliminates the effects of higher harmonics. 

Phase is interpreted in the context of a single waveform as being the difference 
between the start of the cycle (i.e., the zero-crossing point) and the point of zero phase for the 
component of SV(t) of frequency co, expressed as a phase offset. Since the phase offset is an 
10 average over the entire waveform, it may be used as the phase offset from the midpoint of the 
cycle. Ideally, with no zero offset and constant amplitude of oscillation, the phase offset 
should be zero every cycle. In practice, however, it shows a high level of variation and 
provides an excellent basis for correcting mass flow to account for dynamic changes in 
amplitude. 

15 The controller then calculates a phase difference (step 2615). Though a number 

of definitions of phase difference are possible, the analysis assumes that the average phase 
and frequency of each sensor signal is representative of the entire waveform. Since these 
frequencies are different for SVi and SV2, the corresponding phases are scaled to the average 
frequency. In addition, the phases are shifted to the same starting point (i.e., the midpoint of 

20 the cycle on SVi). After scaling, they are subtracted to provide the phase difference. 

The controller next determines the rate of change of the amplitude for the cycle n 
(step 2620): 

amp (end of cycle) - amp (start of cycle) 

roc _ amp « 

period of cycle 

= (amp n+ i- am Pn-i)fr e qn- 

This calculation assumes that the amplitude from cycle n+1 is available when calculating the 

25 rate of change of cycle n. This is possible if the corrections are made one cycle after the raw 

amplitude calculations have been made. The advantage of having an accurate estimate of the 

rate of change, and hence good measurement correction, outweighs the delay in the provision 

of the corrected measurements, which, in one implementation, is on the order of 5 
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milliseconds. The most recently generated information is always used for control of the 
conduit (i.e., for generation of the drive signal). 

If desired, a revised estimate of the rate of change can be calculated after 
amplitude correction has been applied (as described below). This results in iteration to 
5 convergence for the best values of amplitude and rate of change. 



b. Frequency Compensation for Feedback and Velocity Effects 
As noted above, the dynamic aspects of the feedback loop introduce time varying 
shifts in the phase due to the small deviations in amplitude about the set-point. This results 
10 in the measured frequency, which is based on zero-crossings, differing from the natural 
frequency of the conduit. If velocity sensors are used, an additional shift in phase occurs. 
This additional shift is also associated with changes in the positional amplitude of the 
conduit. A dynamic analysis can monitor and compensate for these effects. Accordingly, the 
controller uses the calculated rate of amplitude change to correct the frequency estimate (step 
15 2625). 

The position of an oscillating conduit in a feedback loop that is employed to 
maintain the amplitude of oscillation of the conduit constant may be expressed as: 

X = A(t)sin(<o 0 t-Q(t)), 

where 6 (t) is the phase delay caused by the feedback effect. The mechanical Q of the 
20 oscillating conduit is typically on the order of 1000, which implies small deviations in 
amplitude and phase. Under these conditions, 8(t) is given by: 
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2® 0 A(t) 



Since each sensor measures velocity: 

SV(t) = X(i) = i(Osin[co o ^fl(0]+[co o - 0(0^(0 cos[cd ( 
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25 where Y(t) is the phase delay caused by the velocity effect: 
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Since the mechanical Q of the conduit is typically on the order of 1000, and, hence, 
variations in amplitude and phase are small, it is reasonable to assume: 



CO, 



<M(0 



5 This means that the expression for SV(t) may be simplified to: 

SV{t) » co 0 ^(0 cos(oy - 0(0 - 

and for the same reasons, the expression for the velocity offset phase delay may be simplified 



to: 



Y(0 



A{t) 



<M(0 



10 Summing the feedback and velocity effect phase delays gives the total phase delay: 



cp(0 = 0(0 + Y(0 



A(t) t A(t) 



A(t) 



2® 0 A(t) w 0 A(t) 2® 0 A(t) 



and the following expression for SV(t): 

SV(t) « co 0 ^(?;cos[o) 0 r -<p(t)]. 

From this, the actual frequency of oscillation may be distinguished from the 
natural frequency of oscillation. Though the former is observed, the latter is useful for 
density calculations. Over any reasonable length of time, and assuming adequate amplitude 
control, the averages of these two frequencies are the same (because the average rate of 
change of amplitude must be zero). However, for improved instantaneous density 
measurement, it is desirable to compensate the actual frequency of oscillation for dynamic 
effects to obtain the natural frequency. This is particularly useful in dealing with aerated 
fluids for which the instantaneous density can vary rapidly with time. 

The apparent frequency observed for cycle n is delineated by zero crossings 
occurring at the midpoints of cycles n-1 and n+1. The phase delay due to velocity change 
will have an impact on the apparent start and end of the cycle: 



57 



Attorney Docket No.: 02052-079004 



obs_freq n = 



= obs_freq n _j + 



= obs_freq n _, + 



true 


_freq n 
2n 


(<Pn + . • 


"<P.-l) 


true 


_freq„. 




f 


A A } 








K 4n true _froq n A n+l 4n true_freq n ^4 /I . / ; 


1 


(a 






8n 2 






• 



Based on this analysis, a correction can be applied using an integrated error term: 



1 



f \ 



error sum n = error _ sum n-] -7-5- 



8ti 



\ 



, and 



est_freq n = obs_freq n -error_sum n , 
where the value of error_sum at startup (i.e., the value at cycle zero) is; 



error _ sum o = 



1 



8r 



Ao A, 



Though these equations include a constant term having a value of 1 / 871 , actual data has 
indicated that a constant term of 1 / 871 is more appropriate. This discrepancy may be due to 
unmodeled dynamics that may be resolved through further analysis. 

The calculations discussed above assume that the true amplitude of oscillation, A, 
10 is available. However, in practice, only the sensor voltage SV is observed. This sensor 
voltage may be expressed as: 

SV(t) » (o 0 A(t)cos(a> Q t - <p(t)) 
The amplitude, amp_SV(t), of this expression is: 

amp_SV(t)«co 0 ^;. 

15 The rate of change of this amplitude is: 

roc_amp_SV(t) » (o 0 A(t) 

so that the following estimation can be used: 

k(t) _ roc_amp_SV(t) 
~A(t)~ amp_SV(t) 
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c. Application of Feedback and Velocity Effect Frequency Compensation 

Figs. 27A-32B illustrate how application of the procedure 2600 improves the 
estimate of the natural frequency, and hence the process density, for real data from a meter 
having a one inch diameter conduit. Each of the figures shows 10,000 samples, which are 
5 collected in just over 1 minute. 

Figs. 27 A and 27B show amplitude and frequency data from SVj, taken when 
random changes to the amplitude set-point have been applied. Since the conduit is full of 
water and there is no flow, the natural frequency is constant. However, the observed 
frequency varies considerably in response to changes in amplitude. The mean frequency 
10 value is 8 1 .41 Hz, with a standard deviation of 0.057 Hz. 

Figs. 28A and 28B show, respectively, the variation in frequency from the mean 
value, and the correction term generated using the procedure 2600. The gross deviations are 
extremely well matched. However, there is additional variance in frequency which is not 
attributable to amplitude variation. Another important feature illustrated by Fig. 28B is that 
15 the average is close to zero as a result of the proper initialization of the error term, as 
described above. 

Figs. 29A and 92B compare the raw frequency data (Fig. 29A) with the results of 
applying the correction function (Fig. 29B). There has been a negligible shift in the mean 
frequency, while the standard deviation has been reduced by a factor of 4.4. From Fig. 29B, 
20 it is apparent that there is residual structure in the corrected frequency data. It is expected 
that further analysis, based on the change in phase across a cycle and its impact on the 
observed frequency, will yield further noise reductions. 

Figs. 30A and 30B show the corresponding effect on the average frequency, 
which is the mean of the instantaneous sensor voltage frequencies. Since the mean frequency 
25 is used to calculate the density of the process fluid, the noise reduction (here by a factor of 
5.2) will be propagated into the calculation of density. 

Figs. 31 A and 3 IB illustrate the raw and corrected average frequency for a 2" 
diameter conduit subject to a random amplitude set-point. The 2" flowtube exhibits less 
frequency variation that the 1", for both raw and corrected data. The noise reduction factor is 
30 4.0. 
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Figs. 32A and 32B show more typical results with real flow data for the one inch 
flowtube. The random setpoint algorithm has been replaced by the normal constant setpoint. 
As a result, there is less amplitude variation than in the previous examples, which results in a 
smaller noise reduction factor of 1 .5. 



d. Compensation of Phase Measurement for Amplitude Modulation 
Referring again to Fig. 26, the controller next compensates the phase 
measurement to account for amplitude modulation assuming the phase calculation provided 
above (step 2630). The Fourier calculations of phase described above assume that the 
10 amplitude of oscillation is constant throughout the cycle of data on which the calculations 
take place. This section describes a correction which assumes a linear variation in amplitude 
over the cycle of data. 

Ignoring higher harmonics, and assuming that any zero offset has been 
eliminated, the expression for the sensor voltage is given by: 
15 SV(t) « A x (l + X A t)sm(®t) 

where X A is a constant corresponding to the relative change in amplitude with time. As 
discussed above, the integrals Ii and I 2 may be expressed as: 

2co r- 



I x = — f w SV(t)sin((ot)dt, and 

TT JO 



71 

2x 



I 2 =— f u SV(t)cos((ot)dt. 



2co 

71 

Evaluating these integrals results in: 

f « y 

, and 



I x =A t 



r n \ 



2© 



Substituting these expressions into the calculation for amplitude and expanding as a series in 
results in: 



Amp = A l 



( « 1 ^ 

1 + -*,+—^+... 
V co 8co ) 



Assuming X A is small, and ignoring all terms after the first order term, this may 
25 be simplified to: 
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20 



25 



Amp = A ] 



f 71 ^ 



This equals the amplitude of S V(t) at the midpoint of the cycle (t = n I co) . Accordingly, the 

amplitude calculation provides the required result without correction. 

For the phase calculation, it is assumed that the true phase difference and 
frequency are constant, and that there is no voltage offset, which means that the phase value 
should be zero. However, as a result of amplitude modulation, the correction to be applied to 
the raw phase data to compensate for amplitude modulation is: 



Phase = tan 



X 



Assuming that the expression in brackets is small, the inverse tangent function can be 
ignored. 

A more elaborate analysis considers the effects of higher harmonics. Assuming 
that the sensor voltage may be expressed as: 

SV(t) = (1 H-^^t)[Aj sin(art) + A 2 sin(2<ot) + A 3 sin(ort) + A 4 sin(4c;t)] 

such that all harmonic amplitudes increase at the same relative rate over the cycle, then the 
resulting integrals may be expressed as: 



V (o J 



and 



I 2 — X A (l° A i +80^ 2 +45^3 +32A 4 ) 
60co 



for positive cycles, and 



-1 
60co 



M304 -80^ 2 +45^3 -32A 4 ) 



for negative cycles. 

For amplitude, substituting these expressions into the calculations establishes that 
the amplitude calculation is only affected in the second order and higher terms, so that no 
correction is necessary to a first order approximation of the amplitude. For phase, the 
correction term becomes: 
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60 



/ 30AJ + 80A 2 + 45 A 3 + 32 A 4 N 



\ 



Aj (nX A + cj) 



J 



for positive cycles, and 



-1 
60 



V 



30A 1 -80A 2 H-45A 3 -32A 4 



J 



for negative cycles. These correction terms assume the availability of the amplitudes of the 
higher harmonics. While these can be calculated using the usual Fourier technique, it is also 
possible to approximate some or all them using assumed ratios between the harmonics. For 
example, for one implementation of a one inch diameter conduit, typical amplitude ratios are 
A] = 1 .0, A 2 = 0.01, A 3 = 0.005, and A* = 0.001. 



10 e. Application of Amplitude Modulation Compensation to Phase 

Simulations have been carried out using the digital transmitter, including the 
simulation of higher harmonics and amplitude modulation. One example uses f = 80 Hz, 
Ai(t=0) = 0.3, A 2 = 0, A 3 = 0, A4 = 0, X K = le 5 * 48KHz (sampling rate) = 0.47622, which 
corresponds to a high rate of change of amplitude, but with no higher harmonics. Theory 

15 suggests a phase offset of -0.02706 degrees. In simulation over 1000 cycles, the average 
offset is -0.02714 degrees, with a standard deviation of only 2.17e" 6 . The difference between 
simulation and theory (approx 0.3% of the simulation error) is attributable to the model's 
assumption of a linear variation in amplitude over each cycle, while the simulation generates 
an exponential change in amplitude. 

20 A second example includes a second harmonic, and has the parameters f = 80Hz, 

Ai(t=0) = 0.3, A 2 (t=0) = 0.003, A 3 = 0, A4 = 0, X A = -le* 6 * 48KHz (sampling rate) = 
-0.047622. For this example, theory predicts the phase offset to be +2.706e" 3 , +/-2.66% for 
positive or negative cycles. In simulation, the results are 2.714e J +/ -2.66%, which again 
matches well. 

25 Figs. 33A-34B give examples of how this correction improves real flowmeter 

data. Fig. 33 A shows raw phase data from SVi, collected from a 1" diameter conduit, with 
low flow assumed to be reasonably constant. Fig. 33B shows the correction factor calculated 
using the formula described above, while Fig. 33C shows the resulting corrected phase. The 
most apparent feature is that the correction has increased the variance of the phase signal, 
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while still producing an overall reduction in the phase difference (i.e., SV2 - SVj) standard 
deviation by a factor of 1 .26, as shown in Figs. 34A and 34B. The improved performance 
results because this correction improves the correlation between the two phases, leading to 
reduced variation in the phase difference. The technique works equally well in other flow 
5 conditions and on other conduit sizes. 



f . Compensation to Phase Measurement for Velocity Effect 
The phase measurement calculation is also affected by the velocity effect. A 
highly effective and simple correction factor, in radians, is of the form 

10 c v (t k ) = -ASV(t k ), 

n 

where ASV(tk) is the relative rate of change of amplitude and may be expressed as: 

ISyft ) = SV (h«)-SV(h->) 1 

where tk is the completion time for the cycle for which ASV(tk) is being determined, t^+i is 
the completion time for the next cycle, and is the completion time of the previous cycle. 

15 ASV is an estimate of the rate of change of SV, scaled by its absolute value, and is also 
referred to as the proportional rate of change of SV. 

Figs. 35A-35E illustrate this technique. Fig. 35A shows the raw phase data from 
a single sensor (SVi), after having applied the amplitude modulation corrections described 
above. Fig. 35B shows the correction factor in degrees calculated using the equation above, 

20 while Fig. 35C shows the resulting corrected phase. It should be noted that the standard 
deviation of the corrected phase has actually increased relative to the raw data. However, 
when the corresponding calculations take place on the other sensor (SV2), there is an increase 
in the negative correlation (from -0.8 to -0.9) between the phases on the two signals. As a 
consequence, the phase difference calculations based on the raw phase measurements (Fig. 

25 35D) have significantly more noise than the corrected phase measurements (Fig. 35E). 

Comparison of Figs. 35D and 35E shows the benefits of this noise reduction 
technique. It is immediately apparent from visual inspection of Fig. 35E that the process 
variable is decreasing, and that there are significant cycles in the measurement, with the 
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cycles being attributable, perhaps, to a poorly conditioned pump. None of this is discernable 
from the uncorrected phase difference data of Fig. 35D. 

g. Application of Sensor Level Noise Reduction 
5 The combination of phase noise reduction techniques described above results in 

substantial improvements in instantaneous phase difference measurement in a variety of flow 
conditions, as illustrated in Figs. 36A-36L. Each graph shows three phase difference 
measurements calculated simultaneously in real time by the digital Coriolis transmitter 
operating on a one inch conduit. The middle band 3600 shows phase data calculated using 
10 the simple time-difference technique. The outermost band 3605 shows phase data calculated 
using the Fourier-based technique described above. 

It is perhaps surprising that the Fourier technique, which uses far more data, a 
more sophisticated analysis, and much more computational effort, results in a noisier 
calculation. This can be attributed to the sensitivity of the Fourier technique to the dynamic 
15 effects described above. The innermost band of data 3610 shows the same Fourier data after 
the application of the sensor-level noise reduction techniques. As can be seen, substantial 
noise reduction occurs in each case, as indicated by the standard deviation values presented 
on each graph. 

Fig. 36A illustrates measurements with no flow, a full conduit, and no pump 
20 noise. Fig. 36B illustrates measurements with no flow, a full conduit, and the pumps on. 

Fig. 36C illustrates measurements with an empty, wet conduit. Fig. 36D illustrates 

measurements at a low flow rate. Fig. 36E illustrates measurements at a high flow rate. Fig. 

36F illustrates measurements at a high flow rate and an amplitude of oscillation of 0.03 V. 

Fig. 36G illustrates measurements at a low flow rate with low aeration. Fig. 36H illustrates 
25 measurements at a low flow rate with high aeration. Fig. 361 illustrates measurements at a 

high flow rate with low aeration. Fig. 36J illustrates measurements at a high flow rate with 

high aeration. Fig. 36K illustrates measurements for an empty to high flow rate transition. 

Fig. 36L illustrates measurements for a high flow rate to empty transition. 



30 
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2. Flowtube Level Dynamic Modeling 

A dynamic model may be incorporated in two basic stages. In the first stage, the 
model is created using the techniques of system identification. The flowtube is "stimulated" 
to manifest its dynamics, while the true mass flow and density values are kept constant. The 
response of the flowtube is measured and used in generating the dynamic model. In the 
second stage, the model is applied to normal flow data. Predictions of the effects of flowtube 
dynamics are made for both phase and frequency. The predictions then are subtracted from 
the observed data to leave the residual phase and frequency, which should be due to the 
process alone. Each stage is described in more detail below. 



a. System Identification 
System identification begins with a flowtube full of water, with no flow. The 
amplitude of oscillation, which normally is kept constant, is allowed to vary by assigning a 
random setpoint between 0.05 V and 0.3 V, where 0.3 V is the usual value. The resulting 
15 sensor voltages are shown in Fig. 37A, while Figs. 37B and 37C show, respectively, the 
corresponding calculated phase and frequency values. These values are calculated once per 
cycle. Both phase and frequency show a high degree of "structure." Since the phase and 
frequency corresponding to mass flow are constant, this structure is likely to be related to 
flowtube dynamics. Observable variables that will predict this structure when the true phase 
20 and frequency are not known to be constant may be expressed as set forth below. 

First, as noted above, ASV(tk) may be expressed as: 



SV(t k ,)-SV(t k ,) 



25 This expression may be used to determine ASVi and ASV 2 . 

The phase of the flowtube is related to A, which is defined as ASVi - ASV 2 , 
while the frequency is related to A + , which is defined as ASVi + ASV2. These parameters 
are illustrated in Figs. 37D and 37E. Comparing Fig. 37B to Fig. 37D and Fig. 37C to Fig. 
37E shows the striking relationship between A" and phase and between A + and frequency. 
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Some correction for flowtube dynamics may be obtained by subtracting a multiple 
of the appropriate prediction function from the phase and/or the frequency. Improved results 
may be obtained using a model of the form: 

5 y(k) + ai y(k- 1 ) + . . . + a n y(k-n) = b 0 u(k) + b , u(k- 1 ) + . . . + b m u(k-m), 

where y(k) is the output (i.e., phase or frequency) and u is the prediction function (i.e., A" or 

A + ). The technique of system identification suggests values for the orders n and m, and the 

coefficients a, and bj, of what are in effect polynomials in time. The value of y(k) can be 
10 calculated every cycle and subtracted from the observed phase or frequency to get the 

residual process value. 

It is important to appreciate that, even in the absence of dynamic corrections, the 

digital flowmeter offers very good precision over a long period of time. For example, when 

totalizing a batch of 200kg, the device readily achieves a repeatability of less that 0.03%. 
15 The purpose of the dynamic modeling is to improve the dynamic precision. Thus, raw and 

compensated values should have similar mean values, but reductions in "variance" or 

"standard deviation." 

Figs. 38A and 39A show raw and corrected frequency values. The mean values 

are similar, but the standard deviation has been reduced by a factor of 3.25. Though the 
20 gross deviations in frequency have been eliminated, significant "structure" remains in the 

residual noise. This structure appears to be unrelated to the A + function. The model used is 

a simple first order model, where m = n = 1. 

Figs. 38B and 39B show the corresponding phase correction. The mean value is 

minimally affected, while the standard deviation is reduced by a factor of 7.9. The model 
25 orders are n = 2 and m = 10. Some structure appears to remain in the residual noise. It is 

expected that this structure is due to insufficient excitation of the phase dynamics by setpoint 

changes. 

More effective phase identification has been achieved through further simulation 
of flowtube dynamics by continuous striking of the flowtube during data collection (set point 
30 changes are still carried out). Figs. 38C and 39C show the effects of correction under these 
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conditions. As shown, the standard deviation is reduced by a factor of 3 1 . This more 
effective model is used in the following discussions. 

b. A pplication to Flow Data 

5 The real test of an identified model is the improvements it provides for new data. 

At the outset, it is useful to note a number of observations. First, the mean phase, averaged 
over, for example, ten seconds or more, is already quite precise. In the examples shown, 
phase values are plotted at 82 Hz or thereabouts. The reported standard deviation would be 
roughly 1/3 of the values shown when averaged to 10Hz, and 1/9 when averages to 1 Hz. 

10 For reference, on a one inch flow tube, one degree of phase difference corresponds to about 1 
kg/s flow rate. 

The expected benefit of the technique is that of providing a much better dynamic 
response to true process changes, rather than improving average precision. Consequently, in 
the following examples, where the flow is non-zero, small flow step changes are introduced 

15 every ten seconds or so, with the expectation that the corrected phase will show up the step 
changes more clearly. 

Figs. 38D and 39D show the correction applied to a full flowtube with zero flow, 
just after startup. The ring-down effect characteristic of startup is clearly evident in the raw 
data (Fig. 38D), but this is eliminated by the correction (Fig. 39D), leading to a standard 

20 deviation reduction of a factor of 23 over the whole data set. Note that the corrected 

measurement closely resembles white noise, suggesting most flowtube dynamics have been 
captured. 

Figs. 38E and 39E show the resulting correction for a "drained" flowtube. Noise 
is reduced by a factor of 6.5 or so. Note, however, that there appears to be some residual 
25 structure in the noise. 

The effects of the technique on low (Figs. 38F and 39F), medium (Figs. 38G and 
39G), and high (Figs. 38H and 39H) flow rates are also illustrated, each with step changes in 
flow every ten seconds. In each case, the pattern is the same: the corrected average flows 
(Figs, 39F-39H) are identical to the raw average flows (Figs. 38F-38H), but the dynamic 
30 noise is reduced considerably. In Fig. 39H, this leads to the emergence of the step changes 
which previously had been submerged in noise (Fig. 38H). 
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3. Extensions of Dynamic Monitoring and Compensation Techniques 
The previous sections have described a variety of techniques (physical modeling, 
system identification, heuristics) used to monitor and compensate for different aspects of 
5 dynamic behavior (frequency and phase noise caused by amplitude modulation, velocity 
effect, flowtube dynamics at both the sensor and the flowtube level). By natural extension, 
similar techniques well-known to practitioners of control and/or instrumentation, including 
those of artificial intelligence, neural networks, fuzzy logic, and genetic algorithms, as well 
as classical modeling and identification methods, may be applied to these and other aspects 

10 of the dynamic performance of the meter. Specifically, these might include monitoring and 
compensation for frequency, amplitude and/or phase variation at the sensor level, as well as 
average frequency and phase difference at the flowtube level, as these variations occur within 
each measurement interval, as well as the time between measurement intervals (where 
measurement intervals do not overlap). 

15 This technique is unusual in providing both reduced noise and improved dynamic 

response to process measurement changes. As such, the technique promises to be highly 
valuable in the context of flow measurement. 

I. Aeration 

20 The digital flowmeter provides improved performance in the presence of aeration 

in the conduit. Aeration causes energy losses in the conduit that can have a substantial 
negative impact on the measurements produced by a mass flowmeter and can result in 
stalling of the conduit. Experiments have shown that the digital flowmeter has substantially 
improved performance in the presence of aeration relative to traditional, analog flowmeters. 

25 This performance improvement may stem from the meter's ability to provide a very wide 
gain range, to employ negative feedback, to calculate measurements precisely at very low 
amplitude levels, and to compensate for dynamic effects such as rate of change of amplitude 
and flowtube dynamics, and also may stem from the meter's use of a precise digital 
amplitude control algorithm. 

30 The digital flowmeter detects the onset of aeration when the required driver gain 

rises simultaneously with a drop in apparent fluid density. The digital flowmeter then may 
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directly respond to detected aeration. In general, the meter monitors the presence of aeration 
by comparing the observed density of the material flowing through the conduit (i.e., the 
density measurement obtained through normal measurement techniques) to the known, 
nonaerated density of the material. The controller determines the level of aeration based on 
5 any difference between the observed and actual densities. The controller then corrects the 
mass flow measurement accordingly. 

The controller determines the non-aerated density of the material by monitoring 
the density over time periods in which aeration is not present (i.e., periods in which the 
density has a stable value). Alternatively, a control system to which the controller is 

10 connected may provide the non-aerated density as an initialization parameter. 

In one implementation, the controller uses three corrections to account for the 
effects of aeration: bubble effect correction, damping effect correction, and sensor imbalance 
correction. Figs. 40A-40H illustrate the effects of the correction procedure. 

Fig. 40A illustrates the error in the phase measurement as the measured density 

15 decreases (i.e., as aeration increases) for different mass flow rates, absent aeration correction. 
As shown, the phase error is negative and has a magnitude that increases with increasing 
aeration. Fig. 40B illustrates that the resulting mass flow error is also negative. It also is 
significant to note that the digital flowmeter operates at high levels of aeration. By contrast, 
as indicated by the horizontal bar 4000, traditional analog meters tend to stall in the presence 

20 of low levels of aeration. 

A stall occurs when the flowmeter is unable to provide a sufficiently large driver 
gain to allow high drive current at low amplitudes of oscillation. If the level of damping 
requires a higher driver gain than can be delivered by the flowtube in order to maintain 
oscillation at a certain amplitude, then insufficient drive energy is supplied to the conduit. 

25 This results in a drop in amplitude of oscillation, which in turn leads to even less drive 

energy supplied due to the maximum gain limit. Catastrophic collapse results, and flowtube 
oscillation is not possible until the damping reduces to a level at which the corresponding 
driver gain requirement can be supplied by the flowmeter. 

The bubble effect correction is based on the assumption that the mass flow 

30 decreases as the level of aeration, also referred to as the void fraction, increases. Without 
attempting to predict the actual relationship between void fraction and the bubble effect, this 
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correction assumes, with good theoretical justification, that the effect on the observed mass 
flow will be the same as the effect on the observed density. Since the true fluid density is 
known, the bubble effect correction corrects the mass flow rate by the same proportion. This 
correction is a linear adjustment that is the same for all flow rates. Figs. 40C and 40D 
illustrate, respectively, the residual phase and mass flow errors after correction for the bubble 
effect. As shown, the residual errors are now positive and substantially smaller in magnitude 

than the original errors. 

The damping factor correction accounts for damping of the conduit motion due to 
aeration. In general, the damping factor correction is based on the following relationship 
between the observed phase, <fc,bs, and the actual phase, </w 



Vobs = Wtrue 



1 H ^ 



where 3L is a damping coefficient and k is a constant. Fig. 40E illustrates the damping 
correction for different mass flow rates and different levels of aeration. Fig. 40F illustrates 
the residual phase error after correction for damping. As shown, the phase error is reduced 
substantially relative to the phase error remaining after bubble effect correction. 

The sensor balance correction is based on density differences between different 
ends of the conduit. As shown in Fig. 41, a pressure drop between the inlet and the outlet of 
the conduit results in increasing bubble sizes from the inlet to the outlet. Since material 
flows serially through the two loops of the conduit, the bubbles at the inlet side of the conduit 
(i.e., the side adjacent to the first sensor/driver pair) will be smaller than the bubbles at the 
outlet side of the conduit (i.e., the side adjacent to the second sensor/driver pair). This 
difference in bubble size results in a difference in mass and density between the two ends of 
the conduit. This difference is reflected in the sensor signals <SV, and SV 2 ). Accordingly, 
the sensor balance correction is based on a ratio of the two sensor signals. 

Fig. 40G illustrates the sensor balance correction for different mass flow rates and 
different levels of aeration. Fig. 40H illustrates the residual phase error after applying the 
sensor balance correction. At low flow rates and low levels of aeration, the phase error is 
improved relative to the phase error remaining after damping correction. 

Other correction factors also could be used. For example, the phase angle of each 
sensor signal could be monitored. In general, the average phase angle for a signal should be 
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zero. However, the average phase angle tends to increase with increasing aeration. 
Accordingly, a correction factor could be generated based on the value of the average phase 
angle. Another correction factor could be based on the temperature of the conduit. 

In general, application of the correction factors tends to keep the mass flow errors 
5 at one percent or less. Moreover, these correction factors appear to be applicable over a wide 
range of flows and aeration levels. 



J. Setpoint Adjustment 

The digital flowmeter provides improved control of the setpoint for the amplitude 
10 of oscillation of the conduit. In an analog meter, feedback control is used to maintain the 
amplitude of oscillation of the conduit at a fixed level corresponding to a desired peak sensor 
voltage (e.g., 0.3 V). A stable amplitude of oscillation leads to reduced variance in the 
frequency and phase measurements. 

In general, a large amplitude of oscillation is desirable, since such a large 
15 amplitude provides a large Coriolis signal for measurement purposes. A large amplitude of 
oscillation also results in storage of a higher level of energy in the conduit, which provides 
greater robustness to external vibrations. 

Circumstances may arise in which it is not possible to maintain the large 
amplitude of oscillation due to limitations in the current that can be supplied to the drivers. 
20 For example, in one implementation of an analog transmitter, the current is limited to 100 
mA for safety purposes. This is typically 5-10 times the current needed to maintain the 
desired amplitude of oscillation. However, if the process fluid provides significant additional 
damping (e.g., via two-phase flow), then the optimal amplitude may no longer be sustainable. 
Similarly, a low-power flowmeter, such as the two-wire meter described below, 
25 may have much less power available to drive the conduit. In addition, the power level may 
vary when the conduit is driven by capacitive discharge. 

Referring to Fig. 42, a control procedure 4200 implemented by the controller of 
the digital flowmeter may be used to select the highest sustainable setpoint given a maximum 
available current level. In general, the procedure is performed each time that a desired drive 
30 current output is selected, which typically is once per cycle, or once every half-cycle if 
interleaved cycles are used. 
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The controller starts by setting the setpoint to a default value (e.g., 0.3 V) and 
initializing filtered representations of the sensor voltage (filtered_SV) and the drive current 
(filtered_DC) (step 4205). Each time that the procedure is performed, the controller updates 
the filtered values based on current values for the sensor voltage (SV) and drive current (DC) 
5 (step 4210). For example, the controller may generate a new value for filtered_SV as the 
sum of ninety nine percent of filtered_SV and one percent of SV. 

Next, the controller determines whether the procedure has been paused to provide 
time for prior setpoint adjustments to take effect (step 4215). Pausing of the procedure is 
indicated by a pause cycle count having a value greater than zero. If the procedure is paused, 
10 the controller performs no further actions for the cycle and decrements the pause cycle count 
(step 4220). 

If the procedure has not been paused, the controller determines whether the 
filtered drive current exceeds a threshold level (step 4225). In one implementation, the 
threshold level is ninety five percent of the maximum available current. If the current 
15 exceeds the threshold, the controller reduces the setpoint (step 4230). To allow time for the 
meter to settle after the setpoint change, the controller then implements a pause of the 
procedure by setting the pause cycle count equal to an appropriate value (e.g., 100) (step 
4235). 

If the procedure has not been paused, the controller determines whether the 
20 filtered drive current is less than a threshold level (step 4240) and the setpoint is less than a 
maximum permitted setpoint (step 4245). In one implementation, the threshold level equals 
seventy percent of the maximum available current. If both conditions are met, the controller 
determines a possible new setpoint (step 4250). In one implementation, the controller 
determines the new setpoint as eighty percent of the maximum available current multiplied 
25 by the ratio of filteredjSV to filtered_DC. To avoid small changes in the setpoint (i.e., 
chattering), the controller then determines whether the possible new setpoint exceeds the 
current setpoint by a sufficient amount (step 4255). In one implementation, the possible new 
setpoint must exceed the current setpoint by 0.02 V and by ten percent. 

If the possible new setpoint is sufficiently large, the controller determines if it is 
30 greater than the maximum permitted setpoint (step 4260). If so, the controller sets the 

setpoint equal to the maximum permitted setpoint (step 4265). Otherwise, the controller sets 
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the setpoint equal to the possible new setpoint (step 4270). The controller then implements a 
pause of the procedure by setting the pause cycle count equal to an appropriate value (step 
4235). 

Figs. 43 A-43C illustrate operation of the setpoint adjustment procedure. As 
5 shown in Fig. 43C, the system starts with a setpoint of 0.3 V. At about eight seconds of 

operation, aeration results in a drop in the apparent density of the material in the conduit (Fig, 
43A). Increased damping accompanying the aeration results in an increase in the drive 
current (Fig. 43B) and increased noise in the sensor voltage (Fig. 43C). No changes are 
made at this time, since the meter is able to maintain the desired setpoint. 

10 At about fifteen seconds of operation, aeration increases and the apparent density 

decreases further (Fig. 43 A). At this level of aeration, the driver current (Fig. 43B) reaches a 
maximum value that is insufficient to maintain the 0.3 V setpoint. Accordingly, the sensor 
voltage drops to 0.26 V (Fig. 43C), the voltage level that the maximum driver current is able 
to maintain. In response to this condition, the controller adjusts the setpoint (at about 28 

15 seconds of operation) to a level (0.23 V) that does not require generation of the maximum 
driver current. 

At about 38 seconds of operation, the level of aeration decreases and the apparent 
density increases (Fig. 43 A). This results in a decrease in the drive current (Fig. 43B). At 40 
seconds of operation, the controller responds to this condition by increasing the setpoint (Fig. 
20 43C). The level of aeration decreases and the apparent density increases again at about 48 
seconds of operation, and the controller responds by increasing the setpoint to 0.3 V. 

K. Performance Results 

The digital flowmeter has shown remarkable performance improvements relative 

25 to traditional analog flowmeters. In one experiment, the ability of the two types of meters to 
accurately measure a batch of material was examined. In each case, the batch was fed 
through the appropriate flowmeter and into a tank, where the batch was then weighed. For 
1200 and 2400 pound batches, the analog meter provided an average offset of 500 pounds, 
with a repeatability of 200 pounds. By contrast, the digital meter provided an average offset 

30 of 40 pounds, with a repeatability of two pounds, which clearly is a substantial improvement. 



73 



Attorney Docket No.: 02052-079004 



In each case, the conduit and surrounding pipework were empty at the start of the 
batch. This is important in many batching applications where it is not practical to start the 
batch with the conduit full. The batches were finished with the flowtube full. Some positive 
offset is expected because the flowmeter is measuring the material needed to fill the pipe 
5 before the weighing tank starts to be filled. Delays in starting up, or offsets caused by 
aerated flow or low amplitudes of oscillation, are likely to introduce negative offsets. For 
real batching applications, the most important issue is the repeatability of the measurement. 

The results show that with the analog flowmeter there are large negative offsets 
and repeatability of only 200 pounds. This is attributable to the length of time taken to 
10 startup after the onset of flow (during which no flow is metered), and measurement errors 
until full amplitude of oscillation is achieved. By comparison, the digital flowmeter achieves 
a positive offset, which is attributable to filling up of the empty pipe, and a repeatability is 
two pounds. 

Another experiment compared the general measurement accuracy of the two types 
15 of meters. Fig. 44 illustrates the accuracy and corresponding uncertainty of the 

measurements produced by the two types of meters at different percentages of the meters 1 
maximum recommended flow rate. At high flow rates (i.e., at rates of 25% or more of the 
maximum rate), the analog meter produces measurements that correspond to actual values to 
within 0.15% or less, compared to 0.005% or less for the digital meter. At lower rates, the 
20 offset of the analog meter is on the order of 1.5%, compared to 0.25% for the digital meter. 



L. Self- Validating Meter 

The digital flowmeter may used in a control system that includes self- validating 

sensors. To this end, the digital flowmeter may be implemented as a self-validating meter. 
25 Self-validating meters and other sensors are described in U.S. Patent No. 5,570,300, entitled 

"SELF-VALIDATING SENSORS", which is incorporated by reference. 

In general, a self- validating meter provides, based on all information available to 

the meter, a best estimate of the value of a parameter (e.g., mass flow) being monitored. 

Because the best estimate is based, in part, on nonmeasurement data, the best estimate does 
30 not always conform to the value indicated by the current, possibly faulty, measurement data. 

A self-validating meter also provides information about the uncertainty and reliability of the 
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best estimate, as well as information about the operational status of the sensor. Uncertainty 
information is derived from known uncertainty analyses and is provided even in the absence 
of faults. 

In general, a self-validating meter provides four basic parameters: a validated 
5 measurement value (VMV), a validated uncertainty (VU), an indication (MV status) of the 
status under which the measurement was generated, and a device status. The VMV is the 
meter's best estimate of the value of a measured parameter. The VU and the MV status are 
associated with the VMV. The meter produces a separate VMV, VU and MV status for each 
measurement. The device status indicates the operational status of the meter. 

10 The meter also may provide other information. For example, upon a request from 

a control system, the meter may provide detailed diagnostic information about the status of 
the meter. Also, when a measurement has exceeded, or is about to exceed, a predetermined 
limit, the meter can send an alarm signal to the control system. Different alarm levels can be 
used to indicate the severity with which the measurement has deviated from the 

1 5 predetermined value. 

VMV and VU are numeric values. For example, VMV could be a temperature 
measurement valued at 200 degrees and VU, the uncertainty of VMV, could be 9 degrees. In 
this case, there is a high probability (typically 95%) that the actual temperature being 
measured falls within an envelope around VMV and designated by VU (i.e., from 191 

20 degrees to 209 degrees). 

The controller generates VMV based on underlying data from the sensors. First, 
the controller derives a raw measurement value (RMV) that is based on the signals from the 
sensors. In general, when the controller detects no abnormalities, the controller has nominal 
confidence in the RMV and sets the VMV equal to the RMV. When the controller detects an 

25 abnormality in the sensor, the controller does not set the VMV equal to the RMV. Instead, 
the controller sets the VMV equal to a value that the controller considers to be a better 
estimate than the RMV of the actual parameter. 

The controller generates the VU based on a raw uncertainty signal (RU) that is the 
result of a dynamic uncertainty analysis of the RMV. The controller performs this 

30 uncertainty analysis during each sampling period. Uncertainty analysis, originally described 
in "Describing Uncertainties in Single Sample Experiments," S.J. Kline & F.A. McClintock, 
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Mech. Eng., 75, 3-8 (1953), has been widely applied and has achieved the status of an 
international standard for calibration. Essentially, an uncertainty analysis provides an 
indication of the "quality" of a measurement. Every measurement has an associated error, 
which, of course, is unknown. However, a reasonable limit on that error can often be 
5 expressed by a single uncertainty number (ANSI/ASME PTC 1 9. 1 - 1 985 Part 1 , 
Measurement Uncertainty: Instruments and Apparatus). 

As described by Kline & McClintock, for any observed measurement M, the 
uncertainty in M, w M , can be defined as follows: 

10 where M is true {M true ) with a certain level of confidence (typically 95%). This uncertainty is 
readily expressed in a relative form as a proportion of the measurement (i.e. wm/M). 

In general, the VU has a non-zero value even under ideal conditions (i.e., a 
faultless sensor operating in a controlled, laboratory environment). This is because the 
measurement produced by a sensor is never completely certain and there is always some 

15 potential for error. As with the VMV, when the controller detects no abnormalities, the 

controller sets the VU equal to the RU. When the controller detects a fault that only partially 
affects the reliability of the RMV, the controller typically performs a new uncertainty 
analysis that accounts for effects of the fault and sets the VU equal to the results of this 
analysis. The controller sets the VU to a value based on past performance when the 

20 controller determines that the RMV bears no relation to the actual measured value. 

To ensure that the control system uses the VMV and the VU properly, the MV 
status provides information about how they were calculated. The controller produces the 
VMV and the VU under all conditions-even when the sensors are inoperative. The control 
system needs to know whether VMV and VU are based on "live" or historical data. For 

25 example, if the control system were using VMV and VU in feedback control and the sensors 
were inoperative, the control system would need to know that VMV and VU were based on 
past performance. 

The MV status is based on the expected persistence of any abnormal condition 
and on the confidence of the controller in the RMV. The four primary states for MV status 
30 are generated according to Table 1 . 
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Table 1 



Expected 


Confidence 


MV Status 


Persistence 


inRMV 




not applicable 


nominal 


CLEAR 


not applicable 


reduced 


BLURRED 


short 


zero 


DAZZLED 


long 


zero 


BLIND 



A CLEAR MV status occurs when RMV is within a normal range for given 
5 process conditions. A DAZZLED MV status indicates that RMV is quite abnormal, but the 
abnormality is expected to be of short duration. Typically, the controller sets the MV status 
to DAZZLED when there is a sudden change in the signal from one of the sensors and the 
controller is unable to clearly establish whether this change is due to an as yet undiagnosed 
sensor fault or to an abrupt change in the variable being measured. A BLURRED MV status 

10 indicates that the RMV is abnormal but reasonably related to the parameter being measured. 
For example, the controller may set the MV status to BLURRED when the RMV is a noisy 
signal. A BLIND MV status indicates that the RMV is completely unreliable and that the 
fault is expected to persist. 

Two additional states for the MV status are UNVALID ATED and SECURE. The 

15 MV status is UNVALID ATED when the controller is not performing validation of VMV. 
MV status is SECURE when VMV is generated from redundant measurements in which the 
controller has nominal confidence. 

The device status is a generic, discrete value summarizing the health of the meter. 
It is used primarily by fault detection and maintenance routines of the control system. 

20 Typically, the device status 32 is in one of six states, each of which indicates a different 
operational status for the meter. These states are: GOOD, TESTING, SUSPECT, 
IMPAIRED, BAD, or CRITICAL. A GOOD device status means that the meter is in 
nominal condition. A TESTING device status means that the meter is performing a self 
check, and that this self check may be responsible for any temporary reduction in 

25 measurement quality. A SUSPECT device status means that the meter has produced an 

abnormal response, but the controller has no detailed fault diagnosis. An IMPAIRED device 
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status means that the meter is suffering from a diagnosed fault that has a minor impact on 
performance. A BAD device status means that the meter has seriously malfunctioned and 
maintenance is required. Finally, a CRITICAL device status means that the meter has 
malfunctioned to the extent that the meter may cause (or have caused) a hazard such as a 
5 leak, fire, or explosion. 

Fig. 45 illustrates a procedure 4500 by which the controller of a self-validating 
meter processes digitized sensor signals to generate drive signals and a validated mass flow 
measurement with an accompanying uncertainty and measurement status. Initially, the 
controller collects data from the sensors (step 4505). Using this data, the controller 

10 determines the frequency of the sensor signals (step 4510). If the frequency falls within an 
expected range (step 4515), the controller eliminates zero offset from the sensor signals (step 
4520), and determines the amplitude (step 4525) and phase (step 4530) of the sensor signals. 
The controller uses these calculated values to generate the drive signal (step 4535) and to 
generate a raw mass flow measurement and other measurements (step 4540). 

15 If the frequency does not fall within an expected range (step 45 1 5), then the 

controller implements a stall procedure (step 4545) to determine whether the conduit has 
stalled and to respond accordingly. In the stall procedure, the controller maximizes the driver 
gain and performs a broader search for zero crossings to determine whether the conduit is 
oscillating at all. 

20 If the conduit is not oscillating correctly (i.e., if it is not oscillating, or if it is 

oscillating at an unacceptably high frequency (e.g., at a high harmonic of the resonant 
frequency)) (step 4550), the controller attempts to restart normal oscillation (step 4555) of 
the conduit by, for example, injecting a square wave at the drivers. After attempting to 
restart oscillation, the controller sets the MV status to DAZZLED (step 4560) and generates 

25 null raw measurement values (step 4565). If the conduit is oscillating correctly (step 4550), 
the controller eliminates zero offset (step 4520) and proceeds as discussed above. 

After generating raw measurement values (steps 4540 or 4565), the controller 
performs diagnostics (step 4570) to determine whether the meter is operating correctly (step 
4575). (Note that the controller does not necessarily perform these diagnostics during every 

30 cycle.) 
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Next, the controller performs an uncertainty analysis (step 4580) to generate a raw 
uncertainty value. Using the raw measurements, the results of the diagnostics, and other 
information, the controller generates the VMV, the VU, the MV status, and the device status 
(step 4585). Thereafter, the controller collects a new set of data and repeats the procedure. 
5 The steps of the procedure 4500 may be performed serially or in parallel, and may be 
performed in varying order. 

In another example, when aeration is detected, the mass flow corrections are 
applied as described above, the MV status becomes blurred, and the uncertainty is increased 
to reflect the probable error of the correction technique. For example, for a flowtube 

10 operating at 50% flowrate, under normal operating conditions, the uncertainty might be of 
the order of 0.1 - 0.2% of flowrate. If aeration occurs and is corrected for using the 
techniques described above, the uncertainty might be increased to perhaps 2% of reading. 
Uncertainty values should decrease as understanding of the effects of aeration improves and 
the ability to compensate for aeration gets better. In batch situations, where flow rate 

15 uncertainty is variable (e.g. high at start/end if batching from/to empty, or during temporary 
incidents of aeration or cavitation), the uncertainty of the batch total will reflect the weighted 
significance of the periods of high uncertainty against the rest of the batch with nominal low 
uncertainty. This is a highly useful quality metric in fiscal and other metering applications. 

20 M. Two Wire Flowmeter 

Other embodiments are also contemplated. For example, as shown in Fig. 46, the 
techniques described above may be used to implement a "two-wire" Coriolis flowmeter 4600 
that performs bidirectional communications on a pair of wires 4605. A power circuit 4610 
receives power for operating a digital controller 4615 and for powering the driver(s) 4620 to 

25 vibrate the conduit 4625. For example, the power circuit may include a constant output 

circuit 4630 that provides operating power to the controller and a drive capacitor 4635 that is 
charged using excess power. The power circuit may receive power from the wires 4605 or 
from a second pair of wires. The digital controller receives signals from one or more sensors 
4640. 

30 When the drive capacitor is suitably charged, the controller 461 5 discharges the 

capacitor 4635 to drive the conduit 4625. For example, the controller may drive the conduit 
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once during every 10 cycles. The controller 4615 receives and analyzes signals from the 
sensors 4640 to produce a mass flow measurement that the controller then transmits on the 
wires 4605. 

A number of implementations have been described. Nevertheless, it will be 
5 understood that various modifications may be made. Accordingly, other implementations are 
within the scope of the following claims. 
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